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ABSTRACT 


The  problem  of  determining  the  characteristies  of  the  flow 
of  a  fluid  with  two  free  streamlines  in  a  gravitational  field  is 
investigated.  The  fluid  is  assumed  to  be  incompressible  and  perfect 
and  the  flow  steady  and  irrotational .  Two  cases  are  considered.  The 
first  is  such  that  the  whole  upper  surface  is  a  free  streamline.  The 
second  is  such  that  the  upper  surface  is  partly  covered  by  a  solid 
surface  (S^)  extending  to  infinity  upstream.  This  surface  (S^)  and 
the  one  (S)  on  the  lower  surface  are  prescribed.  The  stream  is 
bounded  downstream  at  a  finite  distance  from  the  point  of  overflow  of 
the  surface  (S)  in  order  to  avoid  a  difficulty  that  occurs  in  the 
analysis  which  occurs  when  the  whole  downstream  flow  is  investigated. 

The  problem  is  reduced  to  the  determination  of  a  complex¬ 
valued  function  in  the  unit  semi-circle  with  non-linear  boundary  conditions. 

The  solution  of  the  problem  is  found  by  the  method  of  successive 
approximations  and  then  proving  that  this  sequence  converges  to  a  unique 
analytic  function  satisfying  the  boundary  conditions.  The  convergence 
is  found  to  depend  upon  the  Froude  number  . 

The  results  are  applied  to  a  case  in  which  the  lower  surface 
(S)  consists  of  a  perfectly  flat  bottom  with  a  vertical  wall. 


A  discussion  of  the  scope  of  the  techniques  follows. 
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CHAPTER  I 


INTRODUCTION 

In  this  thesis  we  will  investigate  a  method  of  solving  a  class 
of  problems  in  two  dimensional  fluid  flow  under  gravity.  In  order  to 
introduce  this  topic,  it  seems  advisable  to  consider  the  equations  with 
which  we  will  be  concerned  as  well  as  some  of  the  main  concepts  of  fluid 
flow.  Then  a  brief  discussion  on  previous  work  will  be  presented. 

In  three  dimensions  the  behavior  of  a  fluid  is  entirely 
determined  by  eleven  equations  in  eleven  unknowns.  These  equations  are 
the  equation  of  continuity,  the  three  equations  of  motion,  the  energy 
equation,  the  equation  of  state,  and  six  equations  relating  the 
components  of  the  stress  tensor  to  the  components  of  the  rate  of 
deformation  tensor.  The  unknown  quantities  are  the  density,  the  three 
components  of  the  velocity,  the  six  components  of  the  stress  tensor, 
and  the  pressure  or  temperature. 

The  mathematical  complexity  is  usually  reduced  by  assuming 
that  (l)  all  stresses  in  the  fluid  are  normal  stresses  and  so  reduce  to 
a  scalar  pressure,  in  which  case  the  fluid  is  said  to  be  ideal  or 
perfect  and  (2)  the  density  is  a  function  of  the  pressure  only,  in 
which  case  the  fluid  is  said  to  be  barotropic. 

The  governing  equations  for  two  dimensional  flow  of  an  ideal 
barotropic  fluid  are  four. 
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(pu)  +  (pv)  = 
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(cont inuity ) 


Su  uSu  vSu  _  .  I 

St  Sx  +  Sy  x  p  Sx  * 


(momentum) 


Sv  uSv  vSv  _  F  .  I  hR 
St  Sx  Sy  y  p  Sy  ' 

p  =  p(p)  , 


(momentum) 


where  p  is  the  density,  u  is  the  x  component  of  the  velocity,  v 
is  the  y  component  of  the  velocity,  F  is  the  x  component  of  body 
force,  Fy  is  the  y  component  of  the  body  force,  p  is  the  pressure, 
and  t  is  the  time. 


The  fluid  motion  is  said  to  be  steady  if  all  of  the  flow  para 
meters  are  independent  of  time. 


Su  Sv  Sp  n 

St  St  St 

The  fluid  motion  is  said  to  be  irrotational  if  the  curl  of  the 
velocity  vector  is  zero;  i.e.,  in  two  dimensions 

Su  Sv  _  n 

Sy  “  Sx  =  u  • 


The  fluid  is  said  to  be  incompressible  if  the  total  derivative 
of  the  density  with  respect  to  the  time  is  zero;  i.e.,  in  two  dimensions 


*£  =  <k  . 

dt  St  Sx 


+  ^£  =  0 
Sy 


When  the  flow  is  irrotational  it  is  well  known  that  there 
exists  a  function  cp  called  the  velocity  potential  which  is  such  that 


. 


In  two  dimensions 
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the  velocity  is  equal  to  the  gradient  of  Cp  . 


If  in  addition  to  being  irrotat ional ,  the  flow  is  incompressible , 


then  the  equation  of  continuity  reduces  to  the  Laplace  equation  for  the 
function  cp  •  Hence  in  two  dimensions,  the  velocity  potential  is  harmonic 
at  all  interior  points  of  the  fluid.  It  is  possible  to  construct  from 
Cp(x,y)  its  harmonic  conjugate  ^/(x,y)  .  This  function  is  called  the 
stream  function  because,  as  it  turns  out,  ifj  is  constant  along  each 


streamline  of  the  flow.  The  function  F(z)  =  Cp(x,y)  +  i^/(x,y)  of  the 
complex  variable  z  =  x  +  iy  is  known  as  the  complex  potential  of  the 
flow.  It  is  regular  at  all  interior  points  of  the  fluid  domain.  As  a 
result  of  the  Cauchy-Riemann  equations  and  the  equations  relating  the 
velocity  components  to  the  velocity  potential 


^cp 

^  ’ 


_  .  M  u  - 
ay  a*  ’  a*  ' 


we  deduce  that 


dF 

W  =  U  -  IV  =  —  , 

dz 


where  w  is  known  as  the  complex  velocity.  Hence  once  F(z)  is  known, 


then  the  velocity  distribution  can  be  found  at  once. 


Usually  the  body  force  F  is  assumed  to  be  conservative  in 


which  case  it  is  well  known  that  there  exists  a  scalar  function  U  such 
that  F*  is  the  gradient  of  (-U)  . 


In  this  work  we  will  assume  that  (l)  the  fluid  is  perfect, 
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barotropic,  and  incompressible,  (2)  the  flow  is  steady  and  irrotat ional , 
(j)  the  body  force  is  conservative,  (4)  the  density  p  is  constant 
and  the  value  of  this  constant  is  1  ,  and  (5)  the  flow  can  be 
restricted  to  a  two  dimensional  flow.  Under  these  assumptions  the  flow 
is  governed  by  the  following  equations, 


(u  v)  _  (Sep  SS) 

'  ’  3x  ’  by  ’ 

sfse  +  is *  0 

2  +  2  ’ 

ax  ay 

TT  P  1/2  2v 

U  +  ■*-  +  — (u  +  v  )  =  constant  . 
p  2 

The  last  equation  is  known  as  the  Bernoulli  equation.  Using 
it  the  pressure  can  be  computed  once  the  velocity  distribution  has  been 
found . 

On  each  boundary  of  the  flow  domain  it  is  necessary  that  some 
condition  be  satisfied  Perhaps  the  pressure  distribution  or  the  normal 
velocity  distribution  (e.g.  a4?/an  =  0  along  a  fixed  solid  boundary) 
is  known.  In  the  case  of  a  streamline  along  which  the  pressure  is 
constant,  the  Bernoulli  equation  supplies  the  boundary  condition.  Such 
a  streamline  is  called  a  free  streamline  and  the  boundary  is  called  a 
free  boundary. 

The  following  is  a  list  of  symbols,  with  their  meanings,  which 
will  be  used  extensively  throughout  the  remainder  of  this  work. 
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p  is  the  density, 
p  is  the  pressure, 

u  is  the  x  component  of  the  velocity, 

v  is  the  y  component  of  the  velocity, 

w  =  u  -  iv  is  the  complex  velocity, 

2  2  1  /2 

V  =  (u  +  v  )  ^  is  the  magnitude  of  the  velocity, 
is  the  velocity  vector, 
is  the  external  body  force, 

U  is  the  gravitational  potential, 

g  is  the  acceleration  due  to  gravity, 
cp  is  the  velocity  potential, 
tp  is  the  stream  function, 

F  is  the  complex  potential, 
w  =  f  +  ig  the  Levy-Civita  function, 

co  =  f  +  ig  the  sequence  of  successive  approximations  of  oo  , 
n  n  n  1  7 

SI  =  F  +  iG  =  co  -  oo,  for  n  >  1  . 
n  n  n  n  1 

Other  symbols  which  may  occur  will  be  defined  as  they  arise. 


CHAPTER  II 


PREVIOUS  WORK 


The  study  of  the  motion  under  gravity  of  a  fluid  with  a  free 
upper  surface  has  drawn  considerable  attention  because  of  its  practical 
significance.  Although  numerical  results  are  available,  the  determination 
of  a  general  mathematical  theory  incorporating  all  of  the  parameters  of 
a  real  flow  is  difficult. 

In  special  cases,  for  example  the  flow  of  a  fluid  over  a  sharp 
angle,  the  flow  of  ideal  fluids  leads  to  theoretical  results  which 
agree  satisfactorily  with  experimental  results. 

H.  Villat  [1]  has  constructed  a  basis  for  the  study  of  the 
motion  under  gravity  of  a  fluid  with  a  free  surface  while  H.  Poncin 
[2]  has  done  a  great  deal  of  work  on  this  problem. 

For  a  review  of  the  known  theorems  to  date  see  Hyers  [5] 
and  Birkhof f  [9] . 

In  his  work  [2]  Poncin  considers  the  steady  irrotational 
flow  of  an  ideal  incompressible  fluid  under  gravity  when  the  flow  is 
supported  by  a  surface  (S)  ,  the  bed  of  the  stream  on  which  is 
specified  the  angle  the  velocity  vector  makes  with  the  horizontal.  The 
flow  also  has  a  single  free  streamline  which  forms  the  upper  surface  of 
the  fluid  domain. 


We  shall  consider  Poncin’s  paper. 


■ 
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In  part  A  the  case  where  the  motion  is  periodic  is  considered,, 
This  does  not  concern  us  directly. 

In  part  B  the  case  where  the  flow  extends  to  infinity  upstream 
and  downstream  is  considered.  The  flow  is  horizontal  and  the  velocity 
is  constant  at  the  extremities  of  the  fluid  domain. 

In  part  C  is  considered  the  generalizations  of  the  problem  in 
part  B,  in  which  the  upper  surface  is  partly  covered  by  a  solid  surface  (S^) 
and  when  on  the  lower  surface  (s),  the  fluid  leaves  the  bed  or  the 
surface  (S)  temporarily,  and  returns  to  it  after  a  short  distance. 

The  description  of  periodic  waves  does  not  concern  us  directly 
in  this  discussion,  but  it  is  interesting  to  note  that  the  problem  reduces 
to  either  solving  an  integral  equation  or  solving  a  mixed  boundary  value 
problem  where  the  boundary  values  are  non-linear.  No  attempt  is  made 
to  solve  the  integral  equation.  The  boundary  value  problem  is  solved 
by  the  method  of  successive  approximations. 

In  part  B  of  the  work,  the  fluid  domain  extends  to  infinity 
upstream  and  downstream  where  the  flow  is  horizontal  and  the  velocity 
assumes  constant  values  V(A)  and  V(B)  respectively.  The  upper 
surface  is  a  free  streamline  and  on  (S),  the  lower  surface,  the  slope 
of  the  velocity  is  prescribed.  The  region  occupied  by  the  fluid  in  the 
complex  potential  plane  ^  is  mapped  conformally  into  the  interior  of 
the  upper  half  of  the  unit  circle  Izl  =  1  of  the  complex  z  plane  by 


mapping 


(2.1) 


. 
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With  the  complex  velocity  given  by  expression 
w(z)  =  V(A)  e 


(2.2) 


where 

w(z)  =  f(x,y)  +  ig(x,y)  , 

the  problem  reduces  to  the  determination  of  an  analytic  function  oj(z) 
which  is  regular  in  the  interior  of  the  upper  unit  semi-circle  and 
satisfies  non-linear  boundary  conditions  on  its  boundary  |z|  =1  .  The 
boundary  value  problem  is  solved  by  the  method  of  successive  approximations. 
(See  the  details  below.)  The  arguments  (x,0)  will  refer  to  the  diameter 
of  the  upper  unit  semi-circle  for  -1  <  x  <  +1  while  the  arguments 
(l,s)  will  refer  to  the  circumference  of  this  unit  semi-circle  when 
r  =  1  and  0  <  s  <  it  for  s  the  arc  length.  On  the  boundary  of  this 
upper  unit  semi-circle  the  function  ljJn(z)  under  the  above  notation  must 
satisfy. 


w  (-1)  =  0  , 

n  ’ 


df 
_ n 

dr 


f^(x,0)  =  @(x)  the  diameter, 

"3g  ,(1,8) 

M-  sin  fn_i(1’s)  e 

(l,s)  =  - — - - - - 

sin  s 


(2.3) 


on  the  circumference  , 


where  |_i  is  a  parameter  equal  to  gQ/(jrV^(A))  and  @(x)  is 
prescribed  on  the  diameter  of  the  unit  semi-circle.  On  the  boundary  of 
this  semi-circle  the  solutions  are 
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r  /-.  \  2  sin  s 

fi(1’s)  -  — — 


g1(l,s)  =  0 


^(x.o)  =  e(x) 


J -1  u  -2u  c 


du 


cos  s+1 


/  1  r ^  e(u)(x  -i)  . 

81  X>  it  J  ^  (x-u) (1-ux) '  U  ’ 

f  (l,s)  =  f  (l,s)  +  &  [  N  (u)  Log  |S!~n  1U+SM^1  du 
nv  lv  it  Jq  n  'sin  (u-sj/21 

8n(l>s)  =  |i  ^  Nn(u)du 
1  J  o 


(2.4) 


where 


fn(x,0)  =  6(x) 


g  (x.o)  =  -  — ^  [  N  (u)  Arctan  f- -  Tan  — }  du 

nv  ’  '  it  /  n  '  ll-x  2J 

J  o 


Sln  f n  1 (1’U)  *3gn  1 (1>“) 

N  (u)  =  - -  e  "*1 

nv  '  sm  u 


In  order  that  oo^(z)  be  well  behaved,  it  is  necessary  that 
0(x)  be  of  the  form 

0(x)  =  @o(x)(l-x2)  (2.5) 

where  Q  (x)  is  bounded  in  -1  <  x  <  1  .  This  means  that  the  flow 
must  be  horizontal  at  infinity  upstream  and  downstream.  The  sequence 
z ) -w  ( z ) }  is  shown  to  be  uniformly  bounded  on  the  boundary  of  the 
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semi-circle  and  uniformly  convergent  to  a  unique  limit  function  which 
is  regular  in  the  interior  of  the  domain  D  and  satisfies  tfye  boundary 
conditions  of  the  problem  providing  p.  satisfies  the  following  relation 


|J.  <  “  exP  (-m+3P  )  ,  3G  =  1-e 
it  mm 


-m 


m  =  2  Arg  Sh  \T3Jj/(2jt)  , 


sin  s 


ds 


(2.6) 


The  pressure  distribution  is  analysed  throughout  the  fluid  domain. 
In  order  that  the  solution  be  acceptable  in  a  hydrodynamical  sense  it  is 
essential  that  the  pressure  be  non-negative  at  all  points  of  the  fluid. 

For  irrotational  flow  of  an  incompressible  ideal  fluid  under  gravity 
with  constant  density,  it  is  evident  from  the  Bernoulli  equation  that 
Ap  <  0  .  This  means  that  the  pressure  is  a  superharmonic  function  in 
the  interior  of  the  fluid  domain.  (See  Courant  and  Hilbert,  "Methods 
of  Mathematical  Physics",  pages  307  and  326.)  By  the  Minimum  Principle 
for  superharmonic  functions,  the  pressure  must  assume  its  minimum  on  the 
boundary  of  this  domain.  It  is  necessary  to  show  that  the  pressure  is 
positive  on  the  boundary  of  the  flow,  in  which  case  it  will  be  positive 
at  all  interior  points. 


For  the  limiting  case  and  for  this  case  only  the  pressure  is 
constant  along  the  free  streamline  (L)  and  this  constant  is  non-negative 
if  the  pressure  infinitely  far  upstream  on  this  streamline  is  non-negative. 
For  any  approximation  °°n(z)  the  difference  between  the  pressure  at  any 
point  on  (l)  and  the  point  infinitely  far  upstream  can  be  made  as 
small  as  possible  providing  n  is  made  large. 


11 


On  (S)  the  study  of  the  pressure  reduces  to  the  study  of  the 
behavior  of  g^(x,0)  .  If  g^(x,0)  is  bounded  then  the  solution  is 
acceptable  in  the  hydrodynamical  sense.  This  is  particularly  evident  if 
@(x)  has  a  discontinuity  of  the  first  kind  at  some  point  of  (S)  .  In 
this  case  g  (x,0)  becomes  infinite  with  the  sign  corresponding  to 
0 (x-O)-0 (x+0)  .  Hence  if  the  lower  surface  (S)  presents  angles  point¬ 
ing  into  the  fluid  region,  the  solution  is  not  acceptable  because  the 
pressure  becomes  negative.  The  case  examined  in  the  work  of  a  surface 
(S)  which  is  horizontal  upstream  and  downstream  but  experiences  a  vertical 
drop  at  some  point  of  the  flow  serves  to  illustrate  this.  At  the  point 
which  experiences  this  vertical  drop,  the  pressure  becomes  negative. 

Further  analysis  shows  that  exterior  to  a  neighbourhood  of  this  point 
whose  radius  is  of  the  order  of  10  cm  ,  the  pressure  is  positive. 

Hence  one  can  find  a  streamline  passing  immediately  outside  this  neighbour¬ 
hood  at  all  points  of  which  the  pressure  is  positive.  In  actual  fact 
the  radius  of  the  above  neighbourhood  is  negligible  in  comparison  to 
the  precision  of  other  measurements  occurring  in  the  problem.  So,  for 
this  reason,  this  streamline  may  be  made  to  represent  the  lower 
boundary  of  the  flow.  The  solution  then  obtained  is  such  that  the 
pressure  is  positive  along  this  boundary  and  hence  inside  the  fluid. 

As  we  shall  see,  if  @(x)  is  continuous  and  different iable* 
then  g  (x,0)  is  well  behaved  and  the  solution  is  valid  in  the 
hydrodynamical  sense. 

In  part  C  a  generalization  is  made-  The  fluid  is  interior  to 
two  horizontal  bands  which  are  a  finite  distance  apart.  Along  these 
contours  the  fluid  moves  on  a  solid  surface  (S),  a  fluid  surface  (G) 


.  * 
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of  the  same  nature  but  at  rest,  or  else  it  is  in  contact  with  the 
atmosphere  (L)  where  the  pressure  is  constant.  The  previous  analysis 
can  easily  be  adapted  to  account  for  cases  arising  in  this  problem. 


In  chapter  one  of  part  C  the  case  of  a  canal  partially  covered 
upstream  by  a  surface  (S^)  is  considered. 


(S) 


Figure  1 

He  deals  specifically  with  the  two  cases  when  Q  as  described  on  this 
upper  surface  is  monotonically  decreasing  and  when  Q  is  such  that  the 
surface  (S^)  is  horizontal  and  experiences  a  vertical  drop  as  indicated 
in  the  diagram  above  (Figure  1). 

The  normal  derivative 

-3g(l,s)  sin  f(l,s) 

^  sin  s 

is  somewhat  more  complicated  than  before  but  poses  no  additional 
difficulties . 

Chapter  three  of  part  C  deals  with  the  case  where  the  contours 
of  the  flow  consist  of  four  arcs.  Refer  to  the  figure  below. 
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(L)  is  the  upper  free  streamline.  (S^)  and  (S^)  form  part  of  the 
solid  flow  bed  (S)  .  On  (S^)  and  (S^)  the  angle  the  velocity  makes 
with  the  horizontal  is  prescribed.  At  B  the  flow  leaves  the  surface 
(s)  and  flows  over  the  fluid  at  rest  between  B  and  C  .  The  dotted 
line  BC  represents  the  streamline.  The  solutions  are  obtained  as  out¬ 
lined  earlier,  the  uniform  boundedness  of  the  sequence  of  successive 
approximations  on  the  boundary  of  the  fluid  as  well  as  its  convergence 
follow  in  a  similar  manner  as  above.  Bounds  on  the  parameter  jj.  are 
obtained  and  the  pressure  distribution  is  obtained  as  before  with  similar 
conclus ions . 

Using  the  results  of  Poncin’ s  work  Dumitrescu,  Ionescu,  and 
Craciun  [10]  adopt  techniques  in  the  study  of  inverse  limits  [14],  [153 
and  [16]  to  study  the  problem  of  determining  the  characteristics  of  the 
flow  bounded  by  a  solid  surface  below  and  a  free  streamline  above  both 
extending  to  infinity  upstream  and  downstream. 

In  general  certain  restrictions  must  be  placed  on  Q  to  obtain 
satisfactory  results.  This  has  been  evident  in  the  papers  discussed. 

Although  the  motion  of  a  fluid  bounded  by  two  free  streamlines 
is  somewhat  more  practical  than  that  bounded  by  one  not  much  work  has 
been  done  on  it  in  comparison  to  the  former.  The  papers  [2],  [14],  [18], 
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[19],  [20],  [21],  [22]  and  [2^]  consider  only  a  single  free  streamline, 

Dumitrescu  and  Ionescu  [ 1 3 ]  have  investigated  the  problem  of 
the  motion  of  an  overflow  under  the  influence  of  gravity.  The  fluid 
is  assumed  to  be  incompressible  and  ideal  and  the  flow  irrotational  and 
steady.  In  order  to  avoid  a  theoretical  difficulty  what  would  occur 
if  the  whole  infinite  flow  were  considered,  the  authors  bound  the 
flow  at  a  finite  distance  downstream  of  the  point  of  overflow.  How¬ 
ever  the  flow  extends  infinitely  far  upstream  where  it  is  horizontal 
with  magnitude  V(A)  .  The  flow  region  (Figure  2)  is  bounded  by  two 
free  streamlines  (L^)  and  (L)  ,  an  equipotential  surface  CC^  and 
a  fluid  bed  (S)  on  which  Q  is  specified. 


The  region  -°o  <  cp  <  ^  0  <  f  <  Q  of  the  complex  potential 

plane  f  =  qp+iy/  representing  the  fluid  domain  is  mapped  into  the  upper 
unit  semi-circle  in  the  z  plane  and  the  boundary  conditions  found. 

The  mapping  and  normal  derivative  of  the  real  part  of  cu(z)  on  the 
circumference  of  the  semi-circle  are 


(1+s)2 

2(l+z2) 


#  I§ir#L  *-38(1’8) 


The  other  boundary  conditions  correspond  to  one  and  two  of  (3°3)°  The 
boundary  condition  above  is  more  complex  than  any  other  one  we  have  seen 
to  this  stage.  Once  the  function  cu(z)  has  been  found  to  be  regular 
in  the  upper  unit  semi-circle  and  to  satisfy  the  boundary  conditions 


then  the  streamlines  can  also  be  found. 


. 
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Utilizing  the  method  of  successive  approximations  the  authors 
solve  the  boundary  value  problem  in  the  same  way  as  Poncin  does*  The 
solution  is  expressed  as  the  limit  of  a  sequence  of  holomorphic  functions 
(^(z)  .  However,  in  examining  the  uniform  boundedness  and  convergence 
of  this  sequence,  the  even  approximations  are  such  that  a  singularity 
always  occurs  at  j(/2  =  s  which  lies  on  the  circumference  of  the  upper 
unit  semi-circle  D  .  The  approximations  are  otherwise  very  well  behaved. 
Due  to  this  fact  one  theoretically  can  deduce  nothing  about  convergence 
of  the  sequence  to  a  solution  of  the  problem.  But  the  theoretical 
results  agree  satisfactorily  with  experimental  results  provided  one 
restricts  himself  to  a  vicinity  of  the  overflow. 

The  problem  is  approached  in  a  slightly  different  manner  in 
a  paper  by  Dumitresu,  lonescu  and  Craciun  [8],  Only  the  semi-infinite 
strip  0  <  ip  <  Q  ,  -°°<cp<0  of  the  complex  potential  plane  is 

mapped  onto  the  new  region.  This  region  is  the  upper  half  of  the  z’ 
plane.  The  surface  (s)  of  the  flow  corresponds  to  the  interval  [l,a] 
where  a  is  determined  by  the  mapping.  The  boundary  conditions  are 
found  and  the  method  of  successive  approximations  is  used  as  above  to 
solve  the  problem.  The  results  will  appear  below  later. 

In  order  to  solve  the  problem  for  approximation  other  than  the 
first,  the  authors  map  the  upper  half  plane  z*  onto  the  upper  unit 
semi-circle  of  the  z  plane.  The  normal  derivative  of  the  function 
F^x'jy')  =  Re  [wn(z 1 ) (z  ’  ) }  is  mapped  into  the  upper  unit  semi¬ 
circle  of  the  z  plane  and  corresponds  to  the  outward  normal  derivative 
on  the  circumference  of  this  semi-circle.  The  solution  is  found  for 
this  semi-circle  and  the  results  are  mapped  back  to  the  z*  plane. 
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I  have  found  the  problem  can  be  simplified  if  one  deals  only 
with  the  semi-circle  and  dispenses  completely  with  the  upper  half  z' 
plane.  The  results  which  I  received  for  the  upper  unit  semi-circle 
were  then  mapped  onto  the  upper  half  plane.  They  are  not  identical  with 
those  obtained  in  reference  [8],  With  the  functions  N(t],z')  and  u)^(z’) 
defined  by  (5.12)  and  (5*3)  the  results  that  I  obtained  for  the  approximation 
to  order  n  were  oo  (z')  equal  to 


Sin  fn-l ( ^ e 


"3gn_i ( °) 


-  \Tt]  (t}-1) 


Log 


N(>z>) 

N(ri,l) 


dr] 


while  those  of  the  authors  of  refernce  [8]  are 

-3g  i  (t),0) 

Sin  fn-l^T1,°^  6  (a-1) 


,  ,,  P  OJ-u  ■Ln_i  \  'b  v  j  e  l  )  Nfn  7  1  'l 

w  (z1  )  +  ^  /  -  ■  - - -  Log  /  ^  J  dr] 

*  -'D  ( 1-  r])  (2T]-a-l)s I  ( T]-l)  (rj-a)  ^ 


The  integrals  in  the  latter  case  do  not  appear  to  be  convergent  while 
those  in  the  former  case  are  convergent  as  we  shall  see  below. 


In  [8]  no  proof  of  the  convergence  of  the  sequence  of  successive 
approximations  is  attempted  although  it  is  stated  explicitly  that  the 
convergence  does  depend  upon  the  parameter  [j.  . 


In  this  paper  the  whole  problem  will  be  re-examined  and  all 
of  the  points  discussed  immediately  above  will  be  clarified.  Although 
an  argument  similar  to  Poncin' s  could  be  used  to  prove  the  uniform 
boundedness  and  uniform  convergence  of  the  sequence  of  successive 
approximations  we  shall  use  a  different  approach  because  the  kernels 
of  the  functions  obtained  here  are  more  complicated  and  require  more 
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work  than  those  used  by  Poncin. 

Some  numerical  computations  are  performed  by  hand  calculators 
and  two  examples  are  partially  examined.  The  first  is  a  vertical  weir 
and  the  second  is  a  type  of  horizontal  jet  emitted  from  an  infinite 
rectangular  box.  The  purpose  of  these  computations  is  to  examine  the 
class  of  flows  to  which  this  theory  is  applicable.  Also  some  results 
are  compared  to  those  of  paper  [8]  in  order  to  determine  the  difference 
in  the  behavior  of  the  results  obtained  in  each  case.  Finally  we  wish 
to  investigate  the  computational  difficulties  if  any  and  determine  the 
feasibility  of  using  an  electronic  computor  to  do  most  of  the 
computation. 


CHAPTER  III 


FORMULATION  AND  OUTLINE  OF  SOLUTION  OF  PROBLEM 


Mathematical  Hypothesis 

Consider  steady  irrotational  flow  in  two  dimensions  of  an 
ideal  incompressible  fluid  under  gravity.  With  reference  to  Figure  2 
the  stream  is  bounded  by  the  streamlines  (L)  and  (L^)  and  the 
equipotential  line  CC*  as  well  as  the  solid  surfaces  (S^)  and  (s) 
which  are  of  course  also  streamlines.  The  upper  surface  (S^)  is 
horizontal  here  and  fixed  but  could  be  arbitrary  apart  from  the  fact 
that  it  becomes  horizontal  at  infinity.  The  lower  surface  is  such  that 
it  tends  toward  the  horizontal  direction  asymptotically  in  the  upstream 
direction. 


Figure  2 


-  19  - 


Let  V(A)  denote  the  velocity  at  infinity  upstream  and  H  the 

fluid  height  at  this  point.  We  will  assume  the  density  is  1  and  that 

the  fluid  remains  in  contact  with  the  surface  (S)  .  If  «  =  i(x  ,y  ) 

o  o 

is  the  equation  of  (S),  then  dy^/dx^  =  slope  V  =  0  on  (s)  . 

After  the  fluid  leaves  point  B,  it  falls  freely  in  a 
gravitational  field  forming  the  free  streamline  (L^)  .  A  similar 
description  applies  to  the  case  of  B*  ,  and  it  forms  the  free  stream¬ 
line  (l)  .  The  stream  is  bounded  at  a  finite  distance  downstream  from 
B  by  the  equipotential  CC*  in  order  to  avoid  a  difficulty  that  would 
occur  in  the  theory  if  the  whole  downstream  flow  were  considered. 

Boundary  Conditions 

Under  the  above  assumptions,  the  problem  is  confined  to  the 

determination  of  a  complex  valued  function  f(x  ,y  )  defined  in  the 

o  o 

fluid  domain  of  figure  2  which  satisfies  the  following  boundary  conditions. 

a)  The  surface  (S)  +  (L^)  is  a  streamline  of  the  flow  and 
we  define  ^/(x^y^)  to  be  zero  there. 

b)  The  surface  (L)  +  (S^)  is  a  streamline  of  the  flow 

so  the  condition  =  Q  must  be  satisfied  there  where  Q  is 

the  quantity  of  fluid  per  second  passing  through  an  equipotential 
surface  one  unit  wide. 

c)  We  will  assume  that  at  a  finite  distance  downstream 
there  exists  an  equipotential  line  along  which  the  velocity  is  constant, 
and  define  Cp(x  ,y  )  there  to  be  zero. 


. 
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d)  The  free  surfaces  (L)  and  (L^)  must  be  determined  in  such 
a  way  that  on  them  the  Bernoulli  equation  is  satisfied.  i0e. 

2 

V  -f  2gyQ  =  constant  .  (3»l) 

The  region  occupied  by  the  flow  is  represented  in  the  complex 
potential  plane  f  =  cp+i^  by  the  lettered  region  ABCC’B’A*  in  Figure  3 
where  the  points  in  Figures  2  and  3  are  made  to  correspond  by  identical 
lettering . 


Figure  5 
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The  transformation 


z*  = 


1  +  exp  jtf/Q 

'V 

1  -  exp  jtf/Q 


F 


(3.2) 


maps  the  region  of  the  complex  potential  f  plane  conformally 

f 

into  the  upper  half  of  the  complex  •>  plane  (Dfl )  .  See  Figures  3 
and  4. 


The  quantities  a  and  b  can  be  approximated  or  they  may  be 
made  arbitrary.  The  points  a  and  3  of  Figure  4  correspond  to  the 
points  -a  and  -b+iQ  of  Figure  3  respectively  and  are  determined 
from  a  and  b  by  the  formulae*. 

2 

f  exp_2a^2_+_l  1 
1  exp  jta/Q  -  1  J  ’ 

(3.3) 

_  f  exp  jtb/Q  -  1  ]£ 

\  exp  jtb/Q  +  1  J 

Note  that  as  b  -»  +oo  ,  3  -»  +1  and  the  whole  upper  surface 

of  the  flow  becomes  a  free  streamline. 

The  region  D  ,  in  the  zr  plane  now  is  mapped  conformally 
z 

onto  the  region  D  consisting  of  the  upper  unit  semi-circle  in  the 

z 

z  plane  in  such  a  way  that  3  and  a  correspond  to  -1  and  +1 
respectively  and  the  strip  [3, a]  corresponds  to  [-1,+1]  .  The 
transformation  is  found  to  bes 

z  =  f  ^  +  i  n (z '  -a)  (z 1  -p )  ^  ^  z»  -  ^  ^  (3.4) 
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where  the  square  root  is  chosen  to  be  real  when  z*  is  real  and  greater 
than  a  .  The  inverse  1:  =  f(z')  of  (3*2)  and  the  inverse  z'  =  z'(z) 
Of  (3.U)  can  be  computed,  and  f(z)  is  found  to  be 


f  (z) 


(CL-P)z  \2 
1+z2  / 


(q-p)z  \2 
z2+i  ; 


+ 


(3-5) 


jX  f—" 

where  we  choose  the  branch  - •—  <  Arg  V  z'  ^  and  th.6  branch  of  th.0 
logarithm  to  be  such  that  -Q  <  Im  f  <  +Q  . 


It  is  convenient  to  introduce  the  Levy-Civita  function  w(z) 

defined  by 


w(z)  =  V(A)  exp  (-ioo(z))  , 


(3.6) 


u(z)  =  f(x,y)  +  lg(x,y)  ,  (3-7) 

where  w(z)  is  the  complex  velocity.  We  must  determine  co(z)  from 
the  boundary  condition  on  the  boundary  of  the  unit  semi-circle  in  the 
z  plane.  Denoting  the  arc  length  on  the  circumference  of  the  semi¬ 
circle  in  the  z  plane  by  s  we  obtain: 


V(l,s)  =  V(A)  exp  g(l,s) 
e(i,s)  =  f(i,s)  , 


(3.8) 


which  when  differentiated  yields: 
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or 


dV(l,s)  =  V(A)g'(l,s)  exp  g(l,s)ds  , 


d  s )  -  *1  exp  (-g(l,s)) 
S  U’S;  "  ds  V(A) 


(3.9) 


From  the  velocity  relationship 


df 

dz 


=  V  exp  (-i@)  , 


(3.10) 


where  ZQ  =  XQ  +  iyQ  denotes  the  original  physical  plane  of  the  flow 
we  obtain 


dz  dx  dy  , 

o  o  ,  o  1  .  „  /  n 

-  =  —  + 1  =  w> exp  > 


df 


dcp 


(3.11) 


along  a  streamline  ip  =  constant. 


Equating  the  imaginary  parts  we  obtain 


V(cp)dy  =  sin  @(cp)dcp 


(3.12) 


Now  differentiate  the  Bernoulli  equation  (3.I) 


VdV  +  gdy  =  0  , 
o 


(3.13) 


and  into  this  substitute  (3*12)  to  obtain 


V^cp)  dV  +  g  sin  0 ( cp ) dcp  =  0 


Substitute  this  into  ( 3 • 9 )  for  gT(l,s)  to  obtain 


.if,  a\  _  _  8  sin  9(<p)  (  (1  u 

8  (1’s'  -  ds  v(a)  ^  }  exp  1  8U’s;; 


(3.14) 


Now  compute  ^  from  equation  (3.5)>  let  Cos  7  =  -(a“P )/(a+(3 ) 
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_  Rp  fd£>  -  ^  dz 


ds 


^  ^dS  ^  Re  X  rlc  ^  > 


dz  ds 


=  Re 


9.fe-g.)  x  _k£ 


rt  /,  2n2 


x 


iz 


(1+z  )  ^  cc+p  +  (g-p)z  a+P  ^  (co-P)z 


1+z 


1+z 


=  Re 


(a-P)  Tan  s 


3t 


{(cc+p 


-2)  Cos  s+a-p)  J 


} 


2  Cos  s 


Vs)  > 


where 


(3.15) 


if  7  >  s  >  3t/2  y 

if  jt  >  s  >  7, 

|  >  s  >  0. 


Substitute  equation  (3*8)  for  @(cp)  and  V(cp)  as  well  as  (3.15)  into 
equations  (3.14)  to  obtain: 


where 


gT(l,s)  =  +p.R1(s)  sin  f(l,s)  exp  (-3g(l,s))  ,  (3  =  16) 


M- 


_ S9_ 

jtV3(A) 


(3=17) 


Fr  being  the  Froude  Number.  Here  R^(s)  and  hence  g?(l,s)  is  zero 
where  s  corresponds  to  the  equipotential  line  CCT  which  is  shown  in 
Figure  2  so  7  >  s  >  jt/2  . 


Now  infinitely  far  upstream  the  velocity  is  V(A)  so  from 
expression  (3*6)  w(z)  at  this  point  must  be  zero.  In  the  z*  plane 


- 
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it  corresponds  to  z*  =  +1  so  in  the  z  plane  it  corresponds  to  z  =  6 
where 

_  g-3-2  n/Ti-3 ) (g-1) 

2-a-p 

for  a  and  (3  are  now  known  quantities. 

Note  as  the  upper  surface  (S^)  recedes  to  infinity,  then 
3=1  and  6  — >  -1  so  in  this  case  the  problem  reduces  to  the  flow 
with  an  upper  free  streamline  extending  to  infinity  upstream. 


(3.18) 


The  problem  can  now  be  formulated  as  follows.  The  function 
w(z)  is  to  be  determined  so  as  to  be  regular  in  the  upper  unit  semi¬ 
circle  of  the  z  plane  and  to  fulfill  the  following  conditions  on  the 
boundary  of  this  semi-circle: 


w(6)  =  0 


5  = 


_  a-3-2  \T ( l-(3 )  (a-1 ) 

2-cc-|3 


f(x,0)  =  @(x )  on  the  diameter  -1  <  x  <  +1 


(3.19) 


g’(l,s)  =  +|iR^(s)  sin  f(l,s)  exp  (-3g(l,s))  on  the  circumference 

0  <  s  <  3t 


where 


R^s)  = 


0 


for  7  >  s  >  — 


2  7 


-  Tan  s 


^a+3  -2  i  \  . 

( - C0S  S+1)  V  ~~n~  +  o 

V.  cc-3  2  2  cos  s 


y  a+p  +  _ad[ 


for  jt>s>7,  ~  >  s  >  0, 


The  first  condition  means  that  at  infinity  upstream  the  velocity 
is  constant  and  is  V(A)  .  The  second  condition  means  that  the  slope 
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of  the  velocity  is  given  on  the  two  surfaces  (s)  and  (S^)  and  is 
represented  on  the  diameter  of  the  upper  unit  semi-circle  by  @(x)  „ 

We  can  assume  that  (S^)  is  perfectly  flat  in  which  case  @(x)  =  0  for 
-1  <  x  <  5  <  0  .  Indeed  we  do  this  in  order  to  make  the  computation 
easier  and  to  lessen  the  difficulties  in  proving  uniform  boundedness 
and  uniform  convergence  of  our  sequence  of  successive  approximations 
which  will  be  obvious  later.  For  the  general  case,  i.e.  for  e(x) 
prescribed  arbitrarly  on  the  surface  (S^)  it  is  easily  seen  in  the 
following  discussion  that  this  reduces  to  the  case  when  3=1  and 
5  =  -1  .  The  third  condition  is  just  the  Bernoulli  equation  for  the 
z  plane. 

Also,  once  oo(z)  is  found  then  the  velocity  distribution 
V(z)  can  be  easily  found  from  equation  (3*6)  and  the  boundaries  of  the 
flow  will  be  determined  by  equation  (3.21)  which  we  now  deduce. 


From  equations  (3*6)  and  (3. 10)  we  obtain 


df 


dzo  =  W)  exp  iw(z)  3?  iz  ’ 


df 


where  ^  is  obtained  from  (3*5).  Hence 


dz  = 


rtV^A) 


I-tV-  x  .  exP  (m(g))„v  x  dz 

(1+z2)2  ( 

1+z 


. (3.20) 


Now  to  obtain  the  streamlines  of  the  flow  and  the  boundaries 
of  the  flow  we  integrate  equation  (3.20).  For  example,  the  boundaries 
of  the  flow  are  found  by  integrating  along  the  boundaries  of  the  unit 


semi -circle . 


* 

nol.»i»,rp9  SJBS-,  j/»k  iron  io 
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ox  2 

4  u^¥ 


ex.P.  (i0J.(x)) _ 

g+P-2  x  \ 
°-P  l+x2/ 


dx 


J-O+P 

2 


„(a-l) 

1+x2 


X 


(3.21) 


Z 

o 


_  — Q— 

itV(A) 


Tan  s  exp  (iu)(l,s)) 


■i  f(^r)  -  + 


When  x^  =  5  and  x  -*  +1  then  the  equation  of  the  surface 
(s)  is  obtained.  When  x^  =  -1  and  x  -* 8  then  the  equation  of  the 
surface  (S^)  is  obtained.  When  s^  =  0  and  s  ->  jt/2  the  equation 
of  the  streamline  (L^)  *-s  obtained.  When  s^  =  jt/2  and  s  -»  7  then 

the  equation  of  the  equipotential  CC’  is  obtained.  When  s^  =  7 
and  s  -*  jt  the  equation  of  the  free  streamline  (L)  is  obtained. 

We  can  use  the  Cauchy  Riemann  equations  for  the  z  plane 

to  write: 


M  (i,s)  =M-  d,s)  =|f  (l,s)  =  g'(l,s)  , 

e 

where  n^  is  the  outward  unit  normal  vector  to  the  region  and  in  this 
case  coincides  with  the  radius. 

The  problem  as  it  is  now  presented  is  such  that  on  one  part 
of  the  boundary  the  real  part  of  w(z)  is  given  and  on  the  other  part 
the  outward  normal  derivative  of  the  real  part  is  given  as  a  function 
of  the  real  and  imaginary  parts  of  w(z)  .  This  problem  has  been  studied 
by  Poncin  and  Villat  [1]  and  [2], 


» 


- 
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The  method  of  successive  approximations  is  used  to  solve  the 
problem.  The  solution  we  seek  will  be  represented  as  the  limit  of  a 
sequence  of  functions  go  (z)  such  that 


go ( z )  =  limit  go  (z)  , 

n  — >  oo 

where  wn(z)  represents  a  function  regular  in  the  upper  unit  semi 
circle  defined  by  the  following  boundary  conditions: 


(0  (5)  =  0  ,  -1  <  S 

n 


a-P-2  \T(l-(3 )  (a-1 ) 

2-a-(3 


f  (x,0)  =  @(x)  on 


df 


~(l,s)  =  +jaR1(s) 


the  diameter  -1  <  x  <  +1  , 

-Jgn_1(l,s) 

sin  f  , (l,s)e  on  the 

n-1 

rt  >  s  >  0 


(3-22) 

c ircumference 


Here  R-^(s)  is  defined  in  (3-19)  ar>d  we  shall  assume  that 

y 

for  I Q  (x)  I  <  t  in  some 
'o'  4 

neighbourhood  of  x  =  6  . 

First  Approximation 

In  solving  this  problem  we  choose  the  first  approximation  to 
be  gravity  free  and  hence  g  =  0  .  Therefore  the  function  go^(z) 
must  be  regular  in  the  unit  semi-circle  and  satisfy: 


e(x)  = 


0  for  -1  <  x  <  5  <  0 


9  (x)  ( x-S )  as  x  -»  6 
o 


+ 


^(6)  =  0 

f^(x,0)  =  @(x)  on  the  diameter  -1  <  x  <  1  , 
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(1,8) 


=  0 


on  the  circumference 


0  <  s  <  it 


The  bounded  solution  to  this  problem  is  known  and  may  be 
found  for  example  in  Poncin  [2]  to  be 


2 

z 


-1 


(x-z) (l-xz) 


(s2-i) 

(x-6) (l-5x) 


dx  . 


(3.23) 


Approximation  to  order  n 

We  construct  the  new  function 

ft  (z)  =  F  (x,y)  +  iG  (x,y)  =  u>  (z)  -  w,(z) 
nv  '  n'  '  n',y'  nv/  lv/ 

which  is  regular  in  the  unit  semi-circle  and  satisfies  the  following 
boundary  conditions. 

Cl  (5)  =  0  -1  <  6  <  0 

nv  —  — 

Fn(x,0)  =0  on  the  diameter  -1  <  x  <  +1  (3.24) 

dF 

(l,s)  =  +jiR1(  s)  sin  fn_1(l,s)  exp  (“3gn.1(l,  s ) ) 
on  the  circumference  0  <  s  <  it 

where  R^(s)  is  defined  by  (3*1 9*)  • 

The  solution  to  this  problem  may  be  found  for  example  in 
H.  Villat  [1]  or  Poncin  [2]  and  it  is 


n  r  *  "3gn_i  ( 1  > s  )  fe1 

^n(z)  =  ^  J  sin  fn_1(l,s)e  \(s)  Log 


is 


-IS  e'l 

-z  e  -o| 
X 


-z 


ls  R, 

e  -o 


ds 


(3.25 
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We  will  investigate  the  behavior  of  F  (x,y)  and  G  (x, y) 

n  n 

on  the  boundary  of  the  unit  semi-circle  but  we  note  that  for  z  =  exp  iu 
we  obtain 


Im 


=  p(u,s)  = 


for 

u  <  s 

for 

u  >  s 

(3.26) 


On  the  boundary  of  the  upper  unit  semi-circle  for  0  <  s  <  it  on  the 
circumference  and  -1  <  x  <  1  on  the  diameter  the  real  and  immaginary 
parts  becomes 


r  / ,  x  2  sin  s  f+1 

fi(1>s)  J  ,  — 

J  -1  u  - 


e(u) 


2u  cos  s+1 


du  , 


/,  ,  i  r+1  e(u)(i-s2)  , 

gl(1’s)  =  «  J_r  ’fu-sjtlfuB)  du  = 


=  constant 


gl(x,°)  =  i  /  ®(u)  {(ui)d-ux)  -  (ui)TI-^B)  }  du  ’ 


(3.27) 


f^XjO)  =  e(x) , 

„  \  u  T*  .  ,  \  ^sn-l^,u^  ,  \  T  i sin  (u+s  )/2 ,  , 

Fn(1’s)  =  *  J  Sln  fn-l(1’u)?  Rl(u)  L0g  Isln  (u-s 5/2-1  du> 


Gnd.s) 


=  ^  /  sin  f  . ( l,u)e 
it  /  n-1'  7 


>it 


-3gn.i(l,u) 


sin  t 


Ri (u) fp(s?u)“2  Arctan  f~ — — 
1 N  v  '  '•cos  u 


r 31  -38n.,(^u)  1+x  u 

G  (x,0)=^  /  sin  f  (l,u)e  2R  (u)  Arctan  (■= -  Tan—)  du  , 

nK  9  '  it  /  n-lv  *  '  lv  '  vl-x  2' 


F  (x,0)  =  0  . 
n 


5} )du  1 
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CHAPTER  IV 

STUDY  OF  SOLUTIONS  OF  THE  PROBLEM 

We  will  study  the  two  cases  3=1  and  3  <  1  with  @(x)  =  0 
for  -1  <  x  <  5  .  In  the  latter  case  the  uniform  convergence  and 
uniform  boundedness  follow  quite  easily  since  the  functions  are  known 
and  well  behaved  in  a  vicinity  of  z  =  -1  .  Otherwise  the  analysis  is 
more  complicated  as  we  shall  see.  The  case  3=1  introduces  a 
singularity  in  the  integrand  which  demands  special  attention.  In  most 
cases  the  lower  surface  (S)  will  become  horizontal  after  some  point 
a  short  distance  upstream,  so  the  fluid  surface  (L)  will  become  horizpntal 
also  at  a  finite  distance  upstream  of  this  point  for  the  case  3=1* 

We  shall  assume  since  the  fluid  is  inviscid  that  the  flow  is  not  affected 
by  imposing  a  flat  surface  (S^)  on  the  upper  surface  of  the  fluid  a 
great  distance  upstream.  This  will  correspond  to  the  case  3  1  and 

will  eleminate  the  necessity  of  treating  a  singularity  in  the  integrand. 

Apart  from  a  constant,  oo^(z)  is  the  same  for  both  cases  so 
we  will  analyse  it  first. 

As  previously  mentioned  we  shall  follow  the  technique  of 
Poncin.  We  shall  however,  postulate  different  conditions  on  Q  and 
use  an  alternative  argument.  As  we  mentioned  in  chapter  two,  Poncin 
considers  the  very  general  case  of  when  Q  is  bounded  and  measurable. 

Our  technique  is  somewhat  restrictive  but  realistic  in  the  physical 
sense.  Then  too,  the  kernels  considered  here  are  quite  different  from 
those  of  Poncin  and  require  careful  investigation. 


'  r  .•  .  uO 
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Study  of  to  (z) 


Clearly  it  suffices  to  study  the  case  0=1.  Then 


fl(x)(z2-l) 

(x-z) (l-xz) 


dx  , 


(4.1) 


is  analytic  outside  the  interval  [ - 1 , +1 ]  of  the  real  axis.  We  shall 
assume  that 


1)  @(x)  is  Holder  continuous  on  [-1,+1]  except  at  a 
finite  number  of  points  where  discontinuities  of  the 
first  kind  occur.  (@  is  Hblder  continuous  in  D* 

if  for  arbitrary  x'  and  x"  in  DT  ,  |@(x' )-0(x") |  < 

A|x’-x"|^  for  A  and  £  <  1  positive  constants.) 

2)  e(x)  =  (i+x)e0(x)  for  |eo(x)  |  <f  ,  (4 

3)  |@(x)-@(l)|  5  ^|x_l|  for  1-e^  -  x  -  ^  ,  where  A  is 

constant  and  is  a  small  number. 


Now  the  equation 


1  r+1 8(x)(ug.i) 

It  7.1  (x-u)(l-ux)  ’ 


C^-3) 


must  be  interpreted  in  the  sense  of  Cauchy  Principal  values.  Clearly, 
however,  if  @(x)  is  Holder  continuous  in  a  vicinity  of  the  point  u 
then  the  Principal  Value  of  (4.3)  exists  (see  [6])  and  is  found  in  terms 
of  the  integral  in  the  ordinary  sense.  However  if  @(x)  presents  any 
jump  discontinuities  then  w(z)  will  have  logarithmic  singularities  at 
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these  points.  For  example  if  (S)  presents  an  angle  pointing  into  the 
stream  or  away  from  the  stream  domain  then  the  velocity  at  this  point 
becomes  infinite  and  zero  respectively. 


Due  to  the  analyticity  of  co  (z)  outside  of  [-1,+1], 
f1(l,s)  is  differentiable  on  (0,jt)  and  if  |@(x)  |  <  2Mjt  then 

|  f  ^ ( 1 ,  s )  |  <  2Hjt  for  0  <  s  <  rt  ,  (4.4) 


because 


2  sin  s 
jf 


dx 


/  \2  2 

(x  -  cos  s )  +  Sin  s 


1 


Now  how 


does  f^(l,s)  behave  as  s  it  and  s  -»  0  ? 


o  '  P  +1 

r-  t-i  \\  ^  2  sin  s  / 

fi(1,s) I  -  ;  J 


(x+1)  0Q(x) 


,  \2  0-2 

(x  -  cos  s )  +  Sin  s 


dx 


<  2M  (-■— S|  Log  (x  -  Cos  s)2  +  (Sin.2  s) 


+1 


-1 


+  ”  (1  4-  Cos  s)) 


=  M  Sin  s  Log  (v  '  ^°S  S )  +  Min:  (l  +  Cos  s) 
L1  +  Cos  sJ 


l-€. 


<  k  (rt-s) 

—  o 


for  constants  k  and  >  0  but  small.  Hence 

o  2 


(^.5) 


f^(l,s)  ->  0  as  s  -*  it, 
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Using  (4.4)  we  can  write 


VM)  -  e(i)| 


le(x)  -  e(D  1 


(x  -  Cos  s)2  +  (Sin  s)2 


dx 


1-e 


2  Sin  s 
Jt 


1  |9(»)-fl(DI 

2  p 

-1  (x  -  Cos  s )  +  (Sin  s) 


dx  + 


-'1-e. 


\eM-e(D\ 


€l  (x  -  Cos  s)d  +  (Sin  s )! 


where  e^  >  0  and  (4.2)  part  4  is  satisfied.  Then 


1-e, 


2  Sin  s 
Jt 


ie(xj-e(i) I 


-1  (x  -  Cos  s)  +  Sin  s 


2 —  dx  <  h  |s  |  , 


for  some  h  >  0  . 


Also 


2  Sin  s 

jt 


|e(x)-e(i) I 


2  2 

1-e^  (x-Cos  s)  +Sin  s 


dx  < 


2A  sin  s 


P1 

m  s  / 

*  A- 


( l-x)dx 


2  2 

e^  (x  -  Cos  s  )  +  Sin  s 


dx 


2A  Sin  s  f  1 


-  g  Log 


2  2 

(l  -  Cos  s )  + Sin  s 


2  2 

(l-e^  -  Cos  s )  +  Sin  s 


PX 

+  (1  -  Cos  s)  —  Arctan  7 - - 7- -rz - ry 

jt  ( e^-2 )  (Cos  s-1) 


e^  Sin  s 


Hence 


f^(l,s)-@(l)  |  -»  0  as  s  -»  0  ,  and 

f^(l,s)  — »  0  ( 1 )  as  s  -*  0 


(b.6) 


dx  j 


. 
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Near  s  =  rc  we  have  the  expression 


1-e, 


f x  ( 1  #  s )  |  <  ko(j(-s), 


while  on  ,  f^(l,s)  is  differentiable  and  continuous  and 

|f^(l,s)  |  <  2Mrt  for  0  <  s  <  jt 

Then  there  exists  a  constant  k^  such  that 

1“«2 

|f1(l,s)  |  <  jjt-s  j  for  0  £  s  <  Jt  » 


(^•7 ) 


Study  of  ^n(z) 


Now  the  real  part  of  the  function  0,  is  zero  on  the  diameter 

n 

of  the  unit  semi-circle  of  the  z  plane  and  is  analytic  in  its 

interior.  The  function  = -ifl^  whose  imaginary  part  is  zero  on  the 
real  axis  and  analytic  inside  the  upper  unit  semi -circle  can  be  continued 
analytically  into  the  lower  unit  ^emi -circle  by  defining  it  there  to  be 
W(z)  .  In  this  way  is  also  extended  to  the  whole  circle  jz[  <  1 

and  Is  analytic  at  all  points  of  this  domain.  Hence  we  will  refer  to 
function  ^n(z)  either  to  the  unit  circle  |z|  =  1  or  just  the  upper 
unit  semi-circle  whenever  it  is  convenient  to  do  so.  Let  D={z||z|<  1}. 

Case  1»  (3  <  1 

The  expression  R  (s)  occurs  in  the  function  Cl  (z)  .  It 


is  defined  by 
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Vs) = 


0 


for  y  >  s  > 


-Tan  s 


sdfci  cos  s+i )  fee.  +  g;P_. .. 

a-(3  J  l  2  2  Cos  s 


—  for  «  >  s  >  7  j 
)  |>s>0. 


The  singularities  at  s  =  y  and  s  =  ^  are  square  root 
singularities  and  hence  are  integrable.  If  @  <  1  and  a  >  1  as  is 
the  case  here  no  other  singularities  occur  in  R  (s)  .  In  the  integrand 
of  ft^(z)  a  singularity  occurs  at  s  =  u  but  this  too  is  integrable 
since  it  is  a  logarithmic  singularity.  It  can  be  shown  that  under  the 
conditions  of  (4.2)  on  @(x)  and  the  boundedness  of  ftR  ^(z)  on  the 
circumference  of  the  circle  Izl  =  1  that  ft  (z)  is  continuous  there 
too.  (See  the  appendix.) 


Using  the  above  facts  fn(l,s)  and  g  (l,s)  will  be  shown 
to  be  uniformly  bounded  on  the  circumference  of  the  circle  |z|  =1  . 
Suppose  that  co^(z)  is  defined  in  the  lower  unit  semi-circle  by  oo^(z) 
where  z  is  in  this  domain.  Since  the  singularities  in  occur  on 

the  diameter  of  the  upper  unit  semi-circle  and  not  on  the  circumference 
f^(l,s)  and  g^(l,s)  are  bounded  on  the  circumference  of  the  circle 
| z |  =  1  .  Let  the  least  upper  bound  be  denoted  by  and 

respectively  for  F.  and  G.  positive  numbers.  Let  F  and  G  be 
upper  bounds  for  the  functions  f  (l,s)  and  gn(l,s)  on  the  circumference 
of  the  circle  |z|  =1  .  The  integrals  in  (3.27)  all  exist  since  the 
singularities  are  of  order  less  than  one  of  the  integrand.  Then 


ro|^ 


' 
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R^u)  Log 


Sin  fu+s)/2 

Sin  (u-s)/2 


du  , 


|glU,s)|  <G  +  ^e 


5G 


n-1 


Sin  u 


|R1(u)  {p(u,s)-2  Arctan  (cos  )  3  1  du 


Using  the  same  argument  as  used  in  proving  wn(z)  continuous 

(see  appendix)  we  can  prove  that  the  integrals  above  are  continuous  on 

the  circumference  of  the  unit  semi-circle.  Here  though  f  ^(l,s)  =  — 

and  g  ,(l,s)  =  0  .  Hence  we  can  take  the  maximum  value  of  these  two 
n-1 

integrals  as  s  traverses  the  circumference  of  the  semi-circle.  Call 

these  maxima  FF  and  GG  respectively  and  define  F  and  G  such 

J  n  n 

that 


5G  , 

F  =  F1  +  ^  e  FF 

n  1  -a 


(4.8a) 


5Gn  1 

G  =  G,  +  ^  e  n"i  GG 
n  1  or 


(4.8b) 


Clearly  G^  >  G^  and  by  the  mathematical  induction  it  follows 
that  G  >  G^  so  the  {G^ }  are  monotonically  increasing  sequences 
of  positive  numbers.  If  these  sequences  are  bounded  above  the  least 
upper  bounds  are  respectively  denoted  by  F  and  G  where  F  and  G 
are  solutions  to  the  equations 

F  =  F.  +  ^  e5G  FF  (4.8c) 

I  it 

G  =  G-  +  ^  e5G  GG 
1  ic 


(4.8d) 
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Suppose  that  we  can  determine  a  value  of  the  parameter  p,  such  that 

equation  (4.8d)  is  satisfied  by  two  values  of  G  which  we  denote  by 

G  and  G  where  G  <  G  .  Since  for  G  >  0 
-L  U  Li  U 

G,  <  G-  +  GG  ^  exp  3G 
1  1  jt  * 

then 


G,  <  G 
1  L  > 


which  is  necessary  since  the  sequence  (G^J  is  monotonically  increasing 


Suppose  for  some  k  ,  G,  is  such  that  G  >  G,  >  G  .  Then  it  is  true 

K  U  K  Lj 


that 


3Gk 

Gk  >  G.  +  *  GG  e  =  G.  . 

K  1  jt  k+1  > 


which  contradicts  the  fact  that  the  sequence  is  increasing.  Hence  G 

La 

is  an  upper  bound  for  the  sequence  [G^J  anc*  sequences  [fn(l,s)} 

and  (g  (l,s)}  are  uniformly  bounded  on  the  boundary  of  the  circle 
Izl  =  1  providing  a  solution  to  (4.8d)  exists,  i.e„, 


^  <  3(GG)  exp  (3G+1) 


(4.9a) 


The  value  of  p.  was  obtained  as  follows.  Let 


yl  =  G, 

y2  =  G1  +  «  GG  6XP 


and  let  us  determine  p.  such  that  these  two  equations  are  satisfied  by 
only  one  value.  This  makes  p.  as  large  as  possible.  Then  y^  must 
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be  tangent  to  the  curve  y^  at  some  point.  Here  the  slope  is  unity  so 

dT  -  1  ■  ^  (GG>  exP  3G 

and  hence 

G  =  3  L°8  3n(GG) 

When  this  value  of  G  is  substituted  in  (4.8d)  an  expression 
similar  to  (4.9a)  is  obtained  but  with  equality. 

The  sequences  {F  (l,s')}  ,  (Gn(l,s))  and  hence  {^n(z)) 

are  uniformly  bounded  on  the  circumference  of  the  circle  |z|  =1  . 

Clearly  the  functions  ft^(z)  are  analytic  in  the  interior  of  this 
circle  and  by  a  theorem  due  to  Montel  [4] 

Theorem  -  "Toute  suite  infinie  de  fonctions  holomorphes 
bornees  dans  un  domaine  (D)  est  generatrice  d’une  suite 
partielle  convergeant  uniformement  vers  une  limite." 

we  know  that  the  sequence  admits  at  least  one  limit  function.  It  remains 

to  find  the  limit  function  and  show  that  it  is  unique  and  that  it  satisfies 

the  boundary  conditions  to  the  problem. 

Now  we  shall  investigate  convergence  of  the  sequence  above. 

We  note  that  since  Ift  -ft  ,  I  =  I  go  -go  I  ,  it  suffices  to  study  the 

1  n  n-1 1  1  n  n-1  1  7 


former.  Note  too  that 


*-i  f  r‘"L 
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-3gn(l,u) 

sin  f  (l,u)e  -  sin  f  . (l,u)e 

n'  7  '  n-1  7  7 


“3g  , (l,u) 

n-i 


"3g  -3g  ,  -3g  , 

<  Isin  f  (e  -e  )  +  e  (sin  f  -  sin  f  )  I 

=  1  n'  '  v  n  n-1  1 


<  e3G  (3|8n-Sn_Ll  +  , 


(4.9b) 


where  G  is  defined  above.  Hence 


|(F3-P2)(l,.)|  <  ±  e*  f*  \h(»)L°S  l§te  [H±||Z||.  (3|G2-Gl|  + 


<  H  e3G  •  4  • 

—  rt 


L  •  FF  , 


where  L  =  Maximum  {  |G  -G  |,  ^2_^1  ^  ^7 


Similarly 


I (G  -G2)(1,s)|  <  ^  e3G  .  v  •  L 


GG 


and  thus  assume  that 


?  -F  I  <  (4  ^  e5GM) 
n  n-1  1  —  jt 


U.  „5GM^n-2  =  en-2  , 


_  I/,  (j.  3Gw\n_2  n-2 

G  -G  i  <  (4  —  e^  M)  =  e„  , 
1  n  n-1  1  —  v  jt  * 


M 


for  all  n  <  N  for  N  given  and  /  is  the  Maximum  of  FF  and  GG  . 
(  i* 

It  follows  easily  that 


and 


IVFN.l  I  <  e*'2 


I  I  N-2 

IVgn-iI  *•* 


x)  ldu 
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Uniform  convergence  of  these  sequences  will  follow  providing 
e*<  1  and  hence 

nc-jjjj  e"3°  •  (4-10) 


So  providing  (4.10)  is  satisfied  by  p,  the  sequences  (Gn(l,s)} 

and  (Fn(l,s)}  are  uniformly  convergent  to  functions  G(l,s)  and  F(l,s) 

which  are  continuous  on  the  circumference  of  the  circle  |z|  =  1  since 

Fn(l,s)  and  Gn(l,s)  are  continuous  there.  Now  let  ft^(z)  and 
2 

ft  (z)  be  two  limit  functions  of  the  sequence  {ftn(z)}  .  Then  these 
two  functions  are  analytic  in  the  domain  D  because  each  fin(z)  is 
analytic  there  and  the  sequence  is  uniformly  convergent.  Then  for 
k  =  1,2  we  can  write 

n 

y  +  ^  > 

i=2 

so  for  n  >  m  in  D  ,  we  have  (n  corresponds  to  k=l  and  m  to  k=2) 


ft  (z)  =  limit 
n  -♦  oo 


n 


|ft1(z)-ft2(z) |  <  limit 


n,m 


00 


y  ini“ni-i 

m+l 


On  the  boundary  of  D  we  can  write 


< 


limit 
n,m  -» oo 


n 


m+l 


1  2 

Since  ft  (z)  -  ft  (z)  is  analytic  throughout  D  and  is  zero  on  the 
boundary  of  D  it  must  be  zero  for  all  z  in  the  domain  D  by  the 


' 
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limit  function  is  unique. 


Are  the  boundary  conditions  of  the  problem  satisfied  by  this 


limit  function.  Clearly  the  conditions  one  and  two  of  (3.24)  are 
satisfied.  From  (3.27)  and  (3.19)  we  can  write 


Due  to  the  uniform  convergence  of  the  sequence  (fn(l,s)}  and  (gn(l,s)} 
on  the  circumference  of  the  circle  Izl  =  1  we  can  write 


-3Sn-l(1’s) 

(sin  fn_1(l,s)e  )R1(s)ds  -i- c 


Limit 

n  ->  00 


g 

=  G(l,s)  =  +  ^  f  sin  f(l,u)e  3§(1>u)r  (u)du  +  c  . 
Jt  /  _  i 


Therefore  on  differentiating  this  expression  we  obtain 


G’(l,s)  =  sin  f(l,s)e  ^g^,S'^R1(s)  , 

jt  1 


which  is  the  third  boundary  condition  of  (3.14)  or  equivalently  (3.25). 

Hence  w(z)  is  the  unique  limit  to  the  sequence  of  successive 
approximations  and  satisfies  the  boundary  conditions. 

Case  2.  p=l. 

With  the  following  notation 


. 
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Rl(u)  = 


0  for  ^  <  u  <  7 


-Tan  u 
1  +  cos  u 


,a  +  1  a  -  1 

(  o  + 


2  cos  u 


)2  for 


0  <  u  <  | 
7  <  u  <  jt 


(4.11) 


,  .  -3gn  ,(l,u) 

®n(u)  =  sin  fn_1(l,u)  e  R^u)  , 


the  real  and  imaginary  parts  of  the  sequence  [ft  (z)}  can  be  written 
in  the  following  manner  on  the  boundary  of  the  upper  unit  semi-circle. 


F  (l,s)  =  ^  f  0  (u)  Log  |Sin  /U+S\/pl  du  y  (4.12) 

n  jt  J  n  'sin  (u-sj/21 


G  (l,s)  =  n  /  $  (u)  du 

n  /  n 

J  Jt 

One  may  proceed  as  indicated  in  the  case  (3  <  1  using  an 
adaption  of  a  technique  developed  by  Poncin  [2].  The  difficulty  here 
is  the  singularity  of  R-^(u)  at  u  =  +jt  .  We  shall  attempt  an 
alternative  proof  below  based  on  a  theorem  by  Golusin  [7]. 

Theorem:  Die  Funktion  f(z)  =  u(r,@)  +  iv(r,@)  sei  in  |z|  <  1 

regular;  ferner  sei  u(r,@)  in  |z|  <  1  stetig  und  genuge  auf 
|z|  =  +1  der  Lipschitz-Bedingung 

|u(e)-u(e' ) |  <  k |e-e'  |a  ,  o  <  a  <  1 

Dann  genugt  f(z)  in  |z|  <  1  einer  komplexen  Lipschitz-Bedingung 

|f(z)-f(z')|  <  K  |z-z'|a 


* 
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Here  the  constant  K  is  dependent  on  the  constants  k  and  a  . 

The  problem  as  we  have  set  it  up  is  to  determine  a  sequence  of 
functions  {^(z)}  which  is  regular  in  the  upper  unit  semi-circle  and 
satisfies  the  following  boundary  conditions. 

fi  (-1)  =  0  , 

n 

dGn( 1, s ) 

— r— -  =  +u,4>  (s)  on  the  circumference  of  the  semi-circle, 

ds  ^  n'  ' 

F  (x,0)  =0  on  the  diameter  of  the  semi-circle, 
n 

We  shall  deal  with  a  different  formulation  of  the  problem 

and  pose  the  question  of  determining  a  sequence  of  functions  W^(z) 

regular  in  the  interior  of  the  circle  Izl  =  1  such  that  W  (-1)  =  0 

ii  n 

and  is  specified  by  the  boundary  conditions 

G  (l,s)  =  p.  P  0  (u)du  for  0  <  s  <  jt  , 

-'jt 

(4.12) 

G  (l,s)  =  G  (1,  2jt  -  s)  for  jt  <  s  <  2jt 
n  n  =  = 

Clearly  W  (z)  coincides  with  -ift  (z)  in  the  upper  semi- 
J  n  n 

circle.  We  must  prove  that  G  ,  the  real  part  of  W  ,  is  continuous 

n  n 

in  | z |  <  1  and  satisfies  a  Lipschitz  condition  on  |z|  =  1  in  order 
to  use  the  theorem  of  Golusin  stated  above.  It  is  not  difficult  to 
show  that  G  satisfies  these  conditions  but  we  shall  illustrate  the 
procedure  by  proving  that  G^  satisfies  a  Lipschitz  condition  on  the 
circumference  Izl  =1  .  Once  this  is  done  it  follows  in  particular  that 


■ 
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ie. 


| f (e  )-f (e  )  |  <  K 


er®2 


.a 


for 


0  <  a  <  1  , 


where  0  <  @^,0^  £  2jt  and  K  =  kcr(a)  0  (See  Golusin). 

We  shall  first  show  that  for  0  <  s  ,Sg  £  J* 

1 

|G2(i, Sf) -G2(l, s2)  |  <  k  |s1-s2 |2  .  (4.14) 

The  extension  to  the  interval  [0,2it]  then  follows  by  symmetry 0  By 
setting  sQ  =  it  and  noting  that  €2(1,3)  =0  at  s  =  it  it  follows 
that  G  (l,s)  is  bounded  and  continuous  for  0  <  s  <  2it  .  From 
previous  considerations  such  as  expression  (4.7)  with  constants  k^ 
and  =  1/2  we  have 

g1(l,s)  =  0  , 

1 

|f1(l,s)  |  <  k  |jt-s  |2 

for  0  <  s  <  it  and  hence  f^(l,s)  is  a  bounded  continuous  function 
on  the  interval  [0,it]  . 

To  prove  (4.14)  we  shall  break  the  discussion  into  two  parts: 

7  <  6^  <  s  <  it  , 

I_  0  <  s  <  5  where  it  >  5  >  5 

2  =  ~  1  11 

With  the  form  of  g^(l,s)  above  it  follows  that  Gg(l,it)  =  0  . 


In  the  region  I  ,  it  follows  that 
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r  2 

jG2(1?s1)-G2(1,s2)  |  =  [x  |  /  sin  fn  1(1,u)R1(u)du  | 

rS2 

<  /  \J~ (st-s)  Rl(s)  jds 

w'  c  . 


<  |ik 


^  L 


ds 


N/Jt-S 


VS1 


<  ji  k 


=  2(^k1h1  |s2-s1 


1 

,2 


where  and 


h-  =  Maximum 

h 


sin  s  jt-s 


1  h-  cos  s  cos  s  ^2 


X  (S^T2 


2^i—  )"2} 

cos  s' 


Now  consider  the  interval  and  denote  the  cosed  intervals 

bounded  by  the  pairs  [0,  jt/2]  ,  [jt/2,  y]  ,  and  [y,  jt]  by  and 

L,  respectively. 

0 


If  s^  and  Sg  are  in  the  result  is  trivial 


integrand  is  null. 


iG2(l,s1)-G2(l,s2)  |  =  0  <  | s1-s2 


1 

,2 


If  s^  and  s2  are  in  then 


since  the 
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where 


|G2(l,s  )-62(L,8  )  I  <  nh2  |  f 

L  “  1  J  S 


du 


si  r  it 

1  ~J  5  '  U 

1 

i  i2 

<  2l^h2  j  s2_sl  I  kl  * 


h  =  Maximum 
Li 


sin  s 


■y*  /(l  • s)  ,0+1 


1  +  cos  s 


/<x+i  a-l\ 

(— —  COS  S  +  —T—J 

cos  s  '  2  2 


If  s^  and  Sg  are  in  the  procedure  is  the  same  and 


lG2^1,Sl' "G2^1,S2^I  <  2M-h5  I  S2~S1  I  kl 


If  s^  is  in  and  s^  is  in  then 


*/2 


|G2(1,s1)-G2(.I,s2)  =  n  |  J  $2(u)du| 


ut  ,2 

<  2fih2  |-  -  sj  kx 


<  2|j.h2  |  s2 -s1  |  kL 


If  s,  and  s.  are  in  L0  and  L_  respectively,  then  as 

12  2J 


above  there  exists  a  constant  h  such  that 

3 


I G2 ( 1 ^ S 1 ) "G2  ^ 1 ’ s2  ^  I  S  2^h3  I  s2 ~s 1  I  kl 


If  s^  and  s2  are  in 


and  respectively,  then  by 


IT 


. 
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a  combination  of  the  arguments  above  there  exists  a  constant  h^  such  that 

1 

Ig2(Mi)-g2(M2)  |  <  2^ki  ls1-s2|2 

Clearly  then,  since  the  regions  1^  and  1^  overlap  we  can 
find  a  number  T  which  depends  only  in  the  integral  above  and  is 
independent  of  the  subscript  n  such  that  for  0  1  si  1  s2  1  * 

1 

|G2(1,sl)-G2(1,s2)|  <  nkiT  |s1-s2 |2  .  (4.16) 

It  is  always  possible  to  determine  T  and  although  it  may  sometimes 

be  difficult.  We  may  now  apply  the  theorem  of  Golusin  and  obtain  the 
result : 

1 

l^2(z1)-fi2(z2)  I  <  ^k1To'(|)  lvz2l2  > 


for 

lZ!L  !  —  1 

and  j  z  |  <1  . 

Also  cr  =  <t 

(f)  is 

independent  of 

the 

subscript 

n  which  in  this 

case  is  2  . 

Since 

G  (l,jt)  =  0  and 

■R 

2 

is  zero  on 

the  diameter  it 

follows  that 

n2(-i) 

=  0  „  In  particular 

^S1  ^S2  2 
|^2  (e  )  -  fig  (e  )|  <  lakjTcr  |  s1~s2  |  , 

for  0  <  s^js2  <  Jt  .  Now  choose  |_i  such  that  for  some  D  >  0 

|a(l+crD)T  exp  3k^  D  v/jt  <  D 
For  such  a  |_i  and  D  it  follows  that 


(J.T  <  D 


' 


. 
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In  order  to  keep  |i  as  large  as  possible  one  would  then 
maximize  the  quotient 


D 


(l-H7D)T  exp  (3^  D\Tjt) 


(lt.1T) 


to  determine  D  .  We  deduce  that 


|F2(1,s1)-F2(1,s2)  |  <  k^D  | s1-s2 


1 

,2 


|G2(l,s1)-G2(l,s2)  |  <  k^D  |sl-s2 


1 

,2 


Now  set  s2  =  rt  and  s^  =  s  so  that  with  our  previous  discussion 
f^(l,s)  and  g^(l,s)  we  can  write  for  0  <  s  <  ^ 


on 


|f2(l,s)  |  <  ki(l-HrD)  (rt-s)  <  k1(l+crD)  \T«  , 


(4.18) 


|g2(l,s)|  <  k^D  (jc-s)  <  k1  D  s/it  , 


so  f  (l,s)  and  g2(l,s)  are  bounded  and  continuous  on  0  <  s  <  ^ 
We  now  repeat  the  argument  above  for  G^(l,s)  „ 


In  the  same  way  as  above  we  show  that 


|G  (l,sL)-G  (l,s2)  |  <  M-k1  ( l+(j.Tcr)  e 


3^k1T  n/tt 


<  k  D  s  -s 
—  i  I  i  2 


1 

,2 


TIS2-S1 


1 

,2 


f2  and  g2  we  deduce  that  G^(l,it)  =  0 


Again  from  the  properties  of 
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and  F^  is  zero  in  the  diameter.  In  the  same  way  as  above  =  0  . 

As  above  from  Golusin 


LS1  lso 

|«5(e  )  '  ^(e  )|  <  k1Dcr  |s]L-s2 


1 

,2 


and 


lG3(M1)-G5(1>s2)|  <  kxD  lvs2 


1 

.2 


|F  (1,s1)-F3(1,s2) |  <  kjD  |s1-s2 


1 

,2 


cr 


In  turn  we  obtdin 


l 

f 1, s )  |  <  k^ (  1+crD)  |jt-s  |2  <  k^(l-KTD)  v/jt 

1 

O  r 

g  (l,s)  |  <  k^D  |jt-s  |  <  k^D  slit 


This  procedure  may  be  repeated  and  hence  for  every  n  and  for  0  S  S1>S2  5  * 
we  have 


S  kiaD  lSrS2 


1 

,2 


lGn^1,Sl^"Gn^1,S2^l  =  kl°  lSl"S2 


1 

,2 


f n(l,  s  )  |  <  k^  1+crD )  nAtt 

gn(i,s)|  <  y/rt 
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It  follows  that  the  sequences  (f^)  arK*  (gn)  as  t'*iey  are 
defined  for  the  circle  |z|  <  1  are  uniformly  bounded  on  the 
circumference  of  this  circle.  Hence  (Wn(z)}  is  uniformly  bounded 
on  the  circumference  of  this  circle  but  W^(z)  Is  analytic  at  all 
points  in  this  circle  so  the  sequence  (Wn(z))  is  uniformly  bounded 
in  the  interior  of  the  circle  |z|  =  1  .  By  a  theorem  due  to  Montel  [4  ] 
which  has  already  been  stated  for  the  case  pel  there  exists  a  sub¬ 
sequence  [W  }  which  is  uniformly  convergent  to  an  analytic  function 
nk 

W(z)  .  Since  W^(z)  coincides  with  -i^n(z)  in  the  upper  unit  semi¬ 
circle  it  follows  that  the  sequence  {^(z)}  also  admits  a  sub¬ 
sequence  which  converges  uniformly  to  an  analytic  function.  It  remains 
to  determine  this  limit  function. 

Note  that  the  boundary  values  Fn(l,s)  and  Gn(l,s)  are 
continuous  functions  on  [0,2jt]  . 

Convergence  of  the  Sequence  {^(z)} 

The  proof  of  the  convergence  of  the  sequences  [Fn(l,s)}  and 
(Gn(l,s)}  and  hence  {^(z)}  on  the  boundary  of  the  domain  D  is  more 
complex  for  the  case  {3  =  1  than  for  the  case  (3  <  1  due  to  the  addition 
of  a  singularity  at  u  =  jt  in  the  function  R^(u)  which  is  defined  in 
(4.11)  and  occurs  in  the  integrand  of  the  expression  for  F^  and  . 

This  integrand  is  also  more  complex  than  those  derived  by  Poncin  in  [2] 
which  we  discussed  earlier  and  there  is  of  course  correspondingly  more 
labor  involved  in  proving  uniform  convergence. 


In  order  to  prove  uniform  convergence  of  the  sequences  involved 


. 
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in  this  analysis  we  shall  resort  once  more  to  the  theorem  of  Golusin. 
This  will  enable  us  to  attack  the  singularity  at  u  =  jt  directly. 

We  shall  first  define  a  number  of  functions. 

Dn(s)  =  Gn(1^s)"Gn_i(1^ s)  for  n  >  2  , 

d2(s)  =  G2^^s)  *  v2(z)  =  -i^2(z)  >  (4.20) 

vn(z)  =  -i{^n(z)-^n_1(z))  for  n  >  2  , 


From  the  expression  (4.16)  we  can  write 

1 

| D2 ( s i ) ” °2 ^ S2 ^  I  =  ^klT  I S2_S1  I2  *  (4.21) 

Since  D  (s)  is  the  real  part  of  the  analytic  function  vn(z)  the 
hypothesis  of  the  theorem  of  Golusin  is  satisfied  where  the  functions 
are  extended  to  the  whole  unit  semi -circle  as  in  (4.13).  This  yields 


1 

2 

| v2 ( 2 1 )  "V2^Z2^  I  5  M-kjTcr  lz1-z2l 

With  the  definitions  F^(l,s)  =  G^(l,s)  =  0  and  since  Fn(l,s)  =  Gn(l,s)  =  0 
for  s  =  jt  this  equation  yields 


1 

p 

Fg(l,s)-F1(l,s)  |  <  iikjTo-  | s -jt  | 

1 

G2(l,s)-G1(ls)  |  <  |akx  |s-7t  |2  T 


(4.22) 


Let  G  denote  the  upper  bound  for  the  sequence  {|g  (l,s)|)  .  From 


' 
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(4.9^)  and  the  expressions  (4.20)  we  have  for 


2 

I ( 8 ! ) ”D3 ( s2 ) I  <  ^  e3G  J'  (3|g2-g1|  +  lf2-fil}  lR1(u)|du  > 

J*l 

s2 

<  H  e5G  J'  {3  |Gg  |  +  P2|)  |Rt(u)  |du  . 


With  (4,22)  and  using  the  value  of  T  as  defined  in  (4.16)  we  obtain 


°2  — 

|D3(s1)-D5(s2)|  <  n  e5G  (a+3)  (iTk1  J*  |R1(u)(«-s)2  |du 


<  fi2T  e3G  (a+3)  kxT  Is^Sgl2 

1 

*  (|iT)2  kL  e3G  (cr+3)  Is^Sg'2 


(4.23) 


The  hypothesis  of  the  theorem  of  Golusin  is  satisfied  and 


hence 


Iv3(z1)-v5(z2)  I  <  ^  e3G  T2a(cr+3)  l^-^l 


1 

,2 


Now  set  z2  = 
we  write 


-1  or  s  =  jf  and  utilizing  the  fact  Gn(l,;t)  =  F^l,jt)  =  0 


G3(1,s)-G2(1,s)|  <  n2  e3G  k^2  (<r+3)(it-«) 

1 

|F3(l,s)-Fg(l,s)  |  <  n2  e3G  kjT2  a(cr+3)  (jt-s )2 


1 

.2 


Proceeding  as  above,  again  we  have  for  s2  > 
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Vsl)"VS2)l  ^  ^  e3G  f  C3|g5-82|  +  |f5-f2|)  |R1(u)|du 


2 

=  U  e3G  J  {3|G3-G2|  +  |F3-F2|)  |R^(u) |du 

S1 

S2  - 

<  \?  e3G  (c T+5)2  e3G  k^2  r  |  (at-s )2  R1(u)|du 

■S 

1 

■zq  2  p 

<  (|a(o-+3)T  )  klM.  |s1-s2|  T  , 


and  hence  from  Golusin 


1 

| ( z i ) ” v4 ( z2 )  I  2  ^(<7+3)2  T5  (e^G)  kj0"  | z1“z2  I2 


so 


1 


|G4(1,s)-G3(1,.)|  <  («-s)£  ^  p2  T, 

1 

|F4(1,s)-F3(1,s)|  <  (*-s)2  |ikL  0-  p2  T 


where 


P  =  n(<r+3)T  e5G  . 


Continuing  in  this  manner  we  obtain  for  all  values  of  the 
subscript  n  >  2, 


1 

Gn(l,s)-Gn_1(l,s)|  <pn'2  |«-s  |2  lik^. 


1 

Fn(1,S)"Fn-l^1,S^  I  -  pn"2  l*'s|2  ^k!Tcr 


(4.24) 


. 
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Hence  providing 


H  < 


1 

T(a+3) 


e 


-3G 


the  sequence  [Wn(z)j  and  hence  the  sequence  {^(z)}  is  uniformly 
convergent  to  an  analytic  function. 


Following  the  same  procedure  as  before  we  know  that  the 
sequence  {^(z)}  converges  to  a  unique  limit  function  which  is  analytic 
in  the  upper  unit  semi-circle  and  satisfies  the  boundary  conditions  of 
the  problem. 
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CHAPTER  V 


EXAMH.ES 


Example  1 

In  the  paper  we  have  referred  to  above  [8],  the  theoretical 
results  are  applied  to  the  case  of  a  perpendicular  weir  with  a  sharp 
corner.  Refer  to  the  figure  below. 


We  shall  refer  to  the  z*  plane  instead  of  the  z  plane  in 
order  to  compare  our  results  with  those  in  reference  [8].  The  surface 
(s)  is  represented  in  the  z’  plane  by  the  segment  [l,a]  .  a  and 
in  this  figure  refer  respectively  to  the  upper  and  lower  edges  of 
the  perpendicular  in  the  zq  plane.  The  function  @(x)  may  thus  be 
written  ^ 


for 

+  V2  for 


1  <  x  <  a 
=  o 

a  <  x  <  a 
o  - 


* 


e(x)  = 


o 


(5.1) 
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a  may  be  determined  from  expresion  3  =  3  and  ot  is  found  to  depend  upon 
the  height  of  the  weir,  The  first  appoximation  of  the  flow  is  given 
by  equation  *(3o23)  for  the  upper  unit  semi -circle,  It  must  be  transformed 
under  expression  (3=4)  to  the  z®  plane.  The  mapping  and  its  inverse 
are  * 


,  a+1  (a-l)z 

z'  "  T  +  ~2  ’ 

1+z 


Z  S 


JL  ,, 

|  +  i{(z'-a)(z«-l)]2  |^2 


t  - 


a+X 


-1 


(5.2) 


1/2 

where  we  shall  choose  the  branch  Arg(z'ra)  '  «  0  for  z'(  real,  z®  >  a 
and  the  branch  Arg  (z'-l)^2  =  0  for  z®  real,  z®  >  1  .  When  (3.23) 
is  transformed  under  this  map  we  obtain ^ 


w1(z')  ■“  {(z’ -l)(z'  -a)}1^2  r 


JSM. 


J  i  {(u-l) (a-u) } 


du 


1/2  z  *  -u 


(5.3) 


As  a  check  it  is  simpler  to  derive  (3=23)  from  (5-3)  since  y 


z  *  -1 ) 

a-uT 


z'  -a 


(u-r 


3 


1  _  (l+x2) (l+z2) 

z *  -u  ~  (z-x5"(l-xz)(a-l) 


du  .  dx 

(x+1  f 


Now  substitute  (5.1)  into  (5=3)  and  integrate  with  a  change  of  variable ^ 


2  a-u 

c  ’ 


. 
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to  obtain 


r  z  *  -co .  1/2 

z  *  -1 


dt 


t2  + 


lzT -1J 


2  _  ^o 

to  a  -1 


(5.4) 


We  shall  choose  the  branches  as  specified  above  and  the  principal  branch 
of  the  logarithm  with  the  result  that  w^(z’)  is  * 


z  *  -a 
z'-l 


1/2 


+  i 


2i  L°8 


a"ao  W2 

a  -1  J 

o  J 


z '  -a 
z'-l 


1/2 


-  i 


a"ao  \l//2 

~  J 

o  J 


(5-5) 


A  branch  point  occurs  at  z1  =  a  but  it  turns  out  that  for  z'  real 

o 

and  l<zf  =  x'  <a. 

=  =  ) 


wl(x' )  =  |  -  |  Log 


a-*'  11/2 


x  *  - 1 


a"ao  l1/2 

a  -1  j 
o  J 


lV2 


a-x 

x»-l 


a"ao  l1/2 

a  -1  J 
o  J 


(5.6) 


The  boundaries  of  the  flow  are  found  from  equation  (5.21)  when 
transformed  to  the  z'  plane.  With  3=1  we  obtain  from  (5.2)^ 


dz'  =  (a-1)  ^~zglp'  dz  j 

(1+z  ) 


so  expression  (3.20)  yields  ^ 


iw(z ' ) 


dz  = 


0  ~  (z'-l) 


dz 1 


Hence  the  boundaries  of  the  flow  for  the  n’th  approximation  are  given  by 
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n 

z 

o 


iw  (x* ) 

n 

e 

n/x*"  (x'-l) 


dx!  , 


(5.7) 


The  approximations  for  n  >  1  are  obtained  from 

for  the.  upper  unit  semi -c ire le  (3.25)  e  Let  s  represent  the  arc  length 
on  the  circumference  of  this  seini-citcle .  We  shall  compute  ^(z*) 
for  the  upper  half  z'  plane  by  transforming  (3.25)  under  the  map  (5.2), 
Define 


q+l  a-1 
2  '  2  Cos  s 


(5.11) 


Then  Cos  s  =  — — ,  Tan  s 

2^-a-l 

specified  the  branches  above. 
(3.25)  becomes  R^r^dr)  where 


.  t!l  where  we  have 

a-1 

Under  the  mapping  (5.2)  R^(s)ds  of 
Rg(T})  is  defined  by  ^ 


/ 


R2(n) 


0 


A  (i-n) 


for  tj  <  0  ^ 

for  0  <  <  1 

a  <  r]  <  00  ^ 


and  the  logarithmic  factor  becomes  Log  N(t),  z’)  where  N  is  defined 
by  the  expression  ^ 


11  II 

(q-1)  (zt-n)^i(2n"a-l)(zt-l)2(zt-a)24-i(2z,-a-l)(n-l)2(ri-a)2  ,  lp 

(q-1 )  (z*  -r])-!-i(2T)-q-l )  (z '  -l)2(z  *  -a)^-i(2z*  -q-1)  ( r)-l)^(  T)-a)^ 

By  our  convention  adopted  above  | N  J  <j:  1  which  agrees  with  the  magnitude 
of  argument  of  the  logarithm  as  it  occurs  in  expressions  (3-25)  and  (3=27)0 
For  example  (z'-l)^-  (z’-q)^^ 


is  positive  for  z1  >  a  and  negative 
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for  z*  <  1  o  The  same  is  true  for  t)  so  Log  jN|  >0  » 

The  expression  is  thus 

n 


0,  (z 

n  ' 


)  = 


J t 


sin  f 


.i(n)e 


-3g 


n 


-1<ri)Log  B(l,)dT) 


(5.13) 


where 


D  =  [0,1]  +  [a,°o)^ 

and 


(5.14) 


If  we  compare  this  expression  with  the  corresponding  expression 
obtained  in  reference  [8]  we  find  that  B(t])  must  be  defined  to  be 


F(t|) 


_ (q-1) _ 

( 1  -  T))  (2  T)-a-l )  T|1,/2  ( rj- 1 ) 1/,£  ( n-a) 1/2 


(5.15) 


We  shall  try  to  find  some  indication  of  how  the  B(t})  as 
defined  in  expressions  (5  =  14)  and  (5. 15)  affect  the  flow.  We  shall 
attempt  several  calculations  by  hand  to  illustrate  the  procedure  in 
evaluating  the  above  integrals  and  to  try  to  achieve  some  insight  into 
their  behavior  and  the  feasibility  of  using  an  automatic  computer. 

The  quantity  a  which  is  the  potential  difference  between 
the  edge  of  the  weir  and  the  equi-poterltial  CC*  is  chosen  to  be 
2Q  .  From  the  expression  (5*3)  a  is  found  to  be  1.0075  •  The  values 
of  the  parameters  H  and  V(A)  are  one  meter  and  ten  meters  per  second, 


i  °*Bori:>  iWM-ihpa  l.,n  jov  9&  \*  *gr,4  9fo 
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so  that  Q  is  ten  cubic  meters  per  second,  the  Froude  number  is 

about  ten  and  the  parameter  ja  is  about  .03  . 


The  height  of  the  weir  is  found  in  conjunction  with  expression 
(5.6)  and  (5.7)  to  be  ^ 


h 


{(x’-l)(a-ao))1//2+{(a-x»)'ao-l)}1'/2  |l/2  dx, 
((x»-l)(a-ao))1//2-{(a-xt  )(aQ-l))1'/2  ^  mQ7(x»-1) 


0 


Here  aQ  appears  in  the  integrand  as  well  as  in  the  limits  of  integration 

and  h  is  known  so  it  is  difficult  to  determine  a  from  h  .  One  must 

o 

substitute  a  sequence  of  aQ  and  compute  the  corresponding  sequence  of 

h  and  then  choose  the  value  of  a  which  closest  approximates  h  . 

o 

In  this  example  <xq  is  chosen  to  be  1,0025  and  the  height  is  computed 
using  the  method  of  tangents  and  the  above  parameters.  The  computation 
is  contained  in  Table  1.  The  height  is  found  to  be  about  ,831  meters 
which  is  less  than  the  height  H  of  the  flow. 


The  remaining  boundaries  of  the  flow  are  found  from  (5.7)  using 
the  expression  for  co^(xf)  as  derived  in  expression  (5.5),  The  result 


(5.8) 


where  P(x’)  is  the  integrand  of  (5*7)  namely  ^ 


P( 


r  v. u. 

x,)  -  {t£ 


(art-l)(x»-a)  ^1/2 


1) (x* -a  )x 

' s  o 


{‘*i 


(a-a  )  (x*  -1)  1/2 


r 

Ik7 


1 ) (x* -a) 


}^r 


(5.9) 


* 


‘  ;  t  ' 
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The  limits  of  integration  x*  and  x'  must  correspond  to  the  stream 
boundary  as  represented  in  the  z'  plane.  For  example,  if  we  wish  to 
find  the  lower  free  streamline  we  set  x’  =  <x  and  vary  x^  along  the 
real  axis  x1  >  a  .  The  expression  (5*8)  is  evaluated  by  the  method 
of  tangents.  In  table  two  the  numerical  computations  are  recorded. 

The  lower  free  streamline  is  found  to  be  higher  than  that  shown  in 
in  figure  6  and  it  becomes  straight  after  one  half  the  distance 
between  B  and  C  .  For  further  information  about  the  boundaries  of 
the  first  approximation  one  should  refer  to  reference  [8],  Appendix  2 
contains  a  brief  discussion  on  the  method  of  tangents. 


Table  1 

Height  of  Weir 


1 

2 

3 

4 

5 

6 

X 

f/x  (x-1)}"1 

A 

dx 

dl 

h 

1.003 

33 2 . 384 

3.731 

.0001 

1.242 

- 

1.004 

249.501 

2.215 

.0001 

.553 

- 

1.005 

199.501 

1.732 

.0001 

.346 

- 

1.006 

166 . I69 

1.447 

.0001 

.241 

- 

1.007 

142.359 

1.211 

.0001 

.172 

.813 

r-A(x-l)1/2  +  (a-x)1/2!1/2 
A=  » 


h  =  -  Z  dl  # 
rt 


Since  the  first  approximation  in  reference  [8]  agrees  with 
those  obtained  here  some  of  the  numerical  results  for  this  approximation 
from  reference  [8]  are  entered  in  tables  one  and  two. 
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Table  2 


Lower  Free  Streamline 


dx 

r2x-2,l/2 

1 

dxl 

dyl 

xi 

yl 

*x-a  J 

\/x  )  (x-1) 

1.01 

.005 

2.828 

99.50 

.053 

.149 

.  .050 

.965 

1.019 

.0125 

1.812 

52.14 

.099 

.182 

.153 

1.144 

1.0375 

.025 

1.581 

26.18 

.116 

.184 

.269 

1.328 

1.0625 

.025 

1.503 

15.52 

.067 

.103 

.337 

1.431 

1.0875 

.025 

1.474 

10.96 

.049 

.072 

.385 

1.503 

1.112 

.025 

1.47 

8.47 

.038 

.056 

.424 

1.559 

1.2 

.15 

1.44 

4.57 

.124 

.178 

.548 

1.737 

1.4 

.25 

1.43 

2,12 

.096 

.138 

.645 

1.875 

1.76 

.47 

1.42 

.99 

.086 

.123 

.731 

1.998 

3. 

2 . 

1.41 

.29 

.107 

.149 

,838 

2.147 

5. 

2. 

1.41 

.11 

.042 

.058 

.879 

2.205 

8. 

4. 

1.41 

.05 

a  037 

.052 

.916 

2.257 

15. 

10. 

1.41 

.02 

.034 

.048 

.950 

2.305 

40. 

40. 

1.41 

.004 

,030 

.042 

.980 

2.347 

100. 

80. 

1.41 

.001 

.015 

.021 

•994 

2.368 

240. 

200. 

1.41 

.0002 

.010 

.0123 

1.005 

2.380 

When  z'  is  real  and  is  not  contained  in  the  interval  [l,a] 
we  can  write  u^(x!)  in  the  form  ^ 


Arc tan 


(a-oo)(*'-l)  -|l/2 
(ao-l)(x'-a)  J  ) 


(5.10) 


which  is  real  if  x*  satisfies  these  conditions.  This  expression  is 
obtained  from  equation  (5.5).  We  can  thus  write 


■  1 

• 

* 

I 

- 

» 

' 


. 
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g1(^S°)  =  o  ^ 


sin  f^(x’,0) 


(a-ao)  (x*  -1) 
(a-1) (x* -a  ) 


o 


When  these  expressions  are  substituted  into  equation  (5.I3)  the  result  is 


) 


«o(l 


)  -if 


f'  f  (a-a0)(n-l)  'll/2 
JD  \  T)(a-X)(ii-%)  J 


Log 


N 


N 


JSiiU 

?rr 


ia_ . 

1-1 


(5.11) 


The  boundaries  of  the  flow  for  the  second  approximation  are 
found  from  expression  (5.7) % 


2 

zo  =  it? 


^  L 


0  IF ,(x')  -0o(x») 


e 


ioo  (x* ) 


Is/x*”  (x*-l) 


dxJ 


(5.12) 


The  quantity  in  the  brackets  in  this  expression  has  been  computed  by  the 

method  of  tangents  for  a  set  of  points  since  it  does  represent  the  flow 

boundaries  for  the  first  approximation,,  Consequently  if  we  evaluate 
2 

zq  at  the  same  set  of  points  then  the  vector  dx^  +  idy^  in  table  2 
representing  the  increments  in  the  flow  boundaries  for  the  first 
approximation  is  essentially  multiplied  by  the  vector  exp  i  «  The 
effect  of  this  operation  is  to  cause  a  rotation  of  the  increment  vector 
dxl  +  idy^  and  an  extension.  On  examination  we  find  the  vector  is 
rotated  an  angle  Fg(x*-)  so  we  can  obtain  some  insight  into  what  is 
going  on  by  examining  the  behavior  for  Fg(x’)  „  One  would  expect  in 
the  case  of  the  free  streamlines  that  Fg(xT)  the  angle  of  rotation  be 
negative  because  gravity  acts  in  such  a  way  as  to  pull  down  the  free 
streamlines  (l)  and  (L^)  from  the  first  approximation  where  gravity 


is  zero. 
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On  examining  (5. 11)  we  find  that  sin  is  always  positive 
in  the  interval  [0,1]  +  [a,°°)  .  Also  by  our  convention  the  logarithmic 
factor  is  always  positive,,  Hence  F^(x')  consists  of  a  positive 
contribution  arising  from  the  interval  [0,1]  and  a  negative 
contribution  arising  from  the  interval  [a,°°)  when  x*  corresponds 
to  the  interval  [a,°°)  .  For  the  lower  free  streamline  we  would 
expect  that  F,_,(x*)  <  ®  • 

If  we  compare  these  results  with  those  of  reference  [8]  we  find 
that  in  using  B(t])  of  (5. 15)  instead  of  (5.14),  the  effect  is  to  scale 
down  the  logarithmic  factor  in  for  i]  >  a  and  increase  it  for 

T)  <  1  .  In  fact  for  large  t]  ,  (5.14)  behaves  like  t]  while 

(5.15)  behaves  like  t)  ^^(a-l)  .  Consequently  the  contribution  from 
the  values  q  >  10  are  significant  for  (5.14)  while  for  (5. 15)  they 
are  not .  So  more  computations  are  necessary  than  those  done  by  reference 
[8]. 


This  discussion  means  that  most  of  the  contribution  to  F^fx’) 
for  x'  >  a  is  coming  from  [0,1]  using  (5.I5)  while  in  the  case  of 
(5.14)  it  is  the  interval  [a,00)  .  Hence  the  results  will  be  such  that 
F^x1  )  will  be  negative  in  one  case  and  positive  in  the  other. 

I  have  computed  F^x’)  For  vari°us  values  of  x!  „  Table  5 
shows  the  results  for  Fg(l5)  »  In  general  when  x*  corresponds  to 
the  upper  free  or  lower  free  streamline  then  F^  results  from  a  negative 
and  a  positive  contribution  as  exhibited  in  table  These  two  contributions 
are  almost  equal  is  magnitude  and  hence  the  sum  is  very  small,,  Consequently 
a  great  number  of  points  should  be  taken  because  accuracy  is  very  important. 
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However  I  did  find  that  was  negative  on  the  interval  [<x,°°)  and 

positive  on  the  interval  [0,1].  When  xl  corresponds  to  the  lower  free 
streamline  then  the  streamline  is  lowering  which  agree  with  the  intuitive 
argument  presented  above.  The  upper  free  streamline  should  be  changed 
also  depending  upon  a  number  of  factors  but  we  have  found  that  for  our 
computation  the  ascent  is  steeper  for  the  second  approximation  than  for 
the  first.  If  one  had  a  large  computor  these  computations  could  be  done 
more  accurately  and  more  quickly  and  the  whole  flow  could  be  analysed. 

But  this  is  not  our  intensions  here.  We  simply  want  to  investigate  this 
method  as  a  means  of  analysing  gravity  flows.  For  more  numerical  results 


see  Appendix  3. 


Table  3 


x  =  15.0 


T) 

At) 

Log  N 

Log  N 

dl 

f (bKti  )(!]-%) 

sTti(t)-1  )(q-ao) 

.02 

.03 

6.19 

.139 

44.57 

1.34 

.05 

.04 

6.17 

.213 

29.36 

1.17 

.1 

.10 

6.11 

.285 

17.91 

1-79  ■ 

.2 

.13 

6.01 

.358 

I6.76 

2.18 

.4 

.2 

5.73 

.380 

15.07 

3.02 

.6 

.2 

5.34 

.311 

17.19 

3.44 

.8 

.2 

4.67 

.180 

25.96 

5.19 

.95 

.09 

3.35 

.050 

67.02 

6.03 

.995 

.01 

1.49 

.006 

243.34 

2.43 

1.01 

.005 

1.10 

-  .0087 

-126.32 

-  .63  . 

1.05 

.0675 

3-33 

-  .050 

-  66.48 

-5.82  . 

1.1 

.04 

3.87 

-  .101 

38.32 

-1.52 

1.55 

.88 

5.71 

-  .683 

-  8.36 

-7.36 

3.0 

2.0 

7. 13 

-  3.462 

-  2.06 

-4.12  . 

6 .0 

4.0 

9.48 

-12.24 

-  .774 

-3.09 

11.0 

4.5 

9.83 

-33.17 

.300 

-1.33 

14.0 

2.5 

11.48 

-48.64 

.246 

-  .59 

16.0 

3.0 

11.63 

-60.00 

.194 

-  .59 

20.0 

17.0 

IO.25 

-84.97 

.121 

-2.06 

50.0 

40.0 

9.25 

-346.5 

.027 

-1.07 

100.0 

100.0 

9.07 

-1000. 

-  .009 

-  .41 

1000. 

5  X  10^ 

8.93 

-31620. 

-2.8  x  10"4 

-I.27 

10000 . 

5  X  104 

8.92 

-1.  x  106 

-8.9  x  10'6 

-  .45 

F0(15)  =  ^  Z  dl  =  -  .037  radians  . 

2  jt 


■ 


Example  2 


We  have  seen  in  the  previous  example  that  it  was  difficult  to 

determine  aQ  for  a  given  value  of  the  weir  height  because  the  height 

was  determined  from  an  integral  equation  where  aQ  appeared  in  the 

limits  of  integration  as  well  as  in  the  integrand.  In  general  if  any 

constants  occur  in  the  definition  of  @(x’)  in  the  interval  [(3, a]  of 

the  complex  z'  plane  or  in  the  interval  [-1,+1]  of  the  complex  z 

plane  then  some  difficulty  of  this  type  will  occur.  The  reason  for 

this  is  that  w  (z’)  in  this  interval  will  be  a  function  of  these 

constants.  The  constants  are  determined  by  the  shape  of  the  surface 

(S)  which  is  known.  We  determine  the  shape  of  (S)  by  the  expression 

z^  for  given  constants.  In  this  expression  the  constants  will  appear 
o 

in  the  integrand  and  in  the  limits  of  integration.  Usually  it  is 
impossible  to  derive  an  expression  for  these  constants  in  terms  of 
known  parameters  in  @(x')  so  some  numerical  method  must  be  used  as 
we  have  found  out  in  the  previous  example. 


As  a  further  illustration  of  the  problems  to  which  this  theory 
is  applicable  and  as  an  illustration  of  how  one  should  interpret  the 
theoretical  results  we  shall  outline  the  approach  to  another  problem 
which  has  interested  many  people. 


Let  us  consider  the  flow  of  a  fluid  through  a  semi-infinite 
pipe  which  is  partially  closed  at  one  end.  See  figure  J . 
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Figure  7 


We  shall  refer  to  the  upper  half  z '  plane  instead  of  the 
unit  semi-circle  because  of  certain  symmetry  conditions  which  will  be 
observed  as  this  discussion  proceeds.  Let  us  define  @(x)  to  be ^ 


e(x)  = 


-*/2 

0 

-Ht/2 


for  (3  <  x  <  3 
=  o 

for  B  <  x  <  a 
=  =0 

for  a  <  x  <  a  . 
o  = 


(5.21) 


The  interval  [3,1]  will  refer  to  the  surface  (S^)  while  the  interval 
[l,a]  will  refer  to  the  surface  (s)  . 


We  shall  use  the  same  numerical  values  for  the  parameters  of 
the  flow  as  was  used  in  the  previous  example.  Hence  the  velocity  V(A) 
is  ten  meters  per  second  and  the  height  H  is  one  meter  so  that  Froude 
number  Fr  is  about  10  and  the  parameter  p.  is  about  .031. 

To  obtain  co^(z’)  we  must  transform  expression  (3.23)  under 
the  mapping  (3-4).  With  Q  defined  above  we  obtain  y 
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( 2 1 ) 


rzT  -a02 

=  l7^j 


,c, 


e(u) 


=  v (z * -a)(z*-py 

^  ^3  */ (u-P)  (a-u) 


da 

z '  -u 


+  iC 


dV 


.2  rz*-a, 

°  •  +fir^ 


(5.22) 

I 

rz‘-PM2  f  dV 

f»’^]  Jo  v2^}  ’ 

Lz 1  -a 


for 


j-  cc-a  ^  1/2 

?i =  j  }  i  and  k 

^  o  J 


PQ-P  -1 1/2 
a-g  J  1 


where  C  is  a  constant  determined  by  the  boundary  condition  ^(1)  =  0 
We  shall  choose  the  branches  Arg  (z * -a) ^ ^  is  zero  for  z*  real  and 
greater  than  a  and  Arg  (z'-p)^^  is  zero  for  z’  real  and  greater 
than  one.  On  taking  the  principal  value  of  the  logarithm  we  obtain 
the  expression  for  co^(z*)^ 


2i 


r  z  *  - 


Log 


T  -a  V 
*-P  J 


+  i 


a -a  -i  2 
o  I 


a  -3 

o 


z 1  -B 
z f  -a 


1 

2 


+  i 


1 

? 


—  -  Log 


z 1  -a 
z 1  -p 


-  i 


a-a  'i  2 

vh 


z T  -3 
z T  -a 


1 

2 


B  -3  1  2 

o 

a-p0  j 


V 


fol  l2 

a-P„  J 


-I-  iC 


The  point  z’  =  a  is  a  branch  point  of  the  first  term  and 
the  point  z’  =  3Q  is  a  branch  point  of  the  second  term.  Consequently 
one  must  be  careful  in  interpreting  this  function  on  the  real  axis. 


It  turns  out  that  four  cases  must  be  considered.  The  appropriate 


results  are 


. 
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Case  1, 


xf  >  a  or 


0  <  x'  <  (3  . 


(xf )  =  Arctan 


-(a-ao)(x,-p)^2 

(a  -P)(x‘-a)j 
-  o 


-  Arctan 


(Po-3)(x'-a)|2 

(a-B  )(**-P>J  + 


Case  2. 


(3  <  x*  <  a  - 

O  0  7 


wi(x* )  =  Jy  Log  lx(x’ )  X  Y(x' )  I  +  iC  ) 


Case  3« 


a  <  x’  <  a 

o 


^  =  2  +  21  L°S  lX^x' )  x  Y(x' ) I  +  iG 


Case  4-. 


p  <  x-  <  Po  , 


co 


1 


|X(xf )  X  Y(x» )  |  +  iC  J 


where 


X(xT ) 


Y(x') 


1 

2 


1 

2 


and  where  C 


is  determined  from  the  boundary  condition 


ux(i) 


(5.24) 


(5.25) 


(5.26) 


(5.27) 


)  (5.28) 


=  0  . 


"•  ■'  .X. 
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C  =  ~  Log  | X ( 1 )  x  Y(l) 


(5.29) 


As  in  the  previous  example  we  shall  choose  a  =  2Q  .  Due  to 
the  conditions  of  symmetry  in  the  surfaces  (s)  and  (S^)  it  turns 


out  that  there  is  also  symmetry  in  the  potentials  of  the  flow  under  no 
gravity.  Hence  we  set  b  =  a  .  a  is  determined  from  expression  (3-5) 
to  be  1.0075  and  3  =  1/a  .  Due  to  the  symmetry  too  we  know  that  the 


potentials  at  the  points  B  and  B*  are  equal.  Then  we  set  =  1.0025 
and  3Q  is  determined  to  be  l/a0  from  (3.3)*  We  must  now  check  to 
see  if  the  heights  h^  and  h^  as  indicated  in  figure  7  total  a 
number  less  than  one. 


Without  the  symmetry  that  we  have  assumed  here  the  constants 


a  ,  a  ,  3  and  3  must  be  altered  in  such  a  way  that  the  computed 
00 

boundaries  coincide  with  the  given  ones.  Sometimes  this  is  difficult 
to  achieve  because  one  does  not  know  in  general  how  any  variation  in 
any  combination  of  these  constants  will  affect  the  boundary  of  the  flow. 
Under  these  symmetry  conditions  we  have  * 


x(l/x')  =  Y(x')  J 


(5.J0) 


SO 


C  =  Log  |X(1) I  =  Log  5.813 


The  formula  for  h^  is  found  using 


as  derived  in  case  3» 


i  f  8  :  * 


J 


ioo  (x? ) 


h  -  =  i  {  -  f  — — - — 

1  "l'  V  J7  (k'-i) 


dx' 


(5.31) 


5.813  it 


|X(x' )  X  Y(x» ) |  I  dx* 


a 


x  *  -1 


Under  the  substitution  x*  =  1/y’  in  this  equation  we  obtain 
an  equation  which  is  idential  with  the  expression  for  -h^  .  This 
means  that  h^  =  -h^  as  one  should  expect.  Using  the  method  of  tangents 
the  height  h^  can  be  computed.  The  computation  is  analogous  to  that 
shown  in  table  1  of  the  first  example.  We  use  the  same  set  of  points 
also.  The  height  is  found  to  be  about  .3  meters  which  is  clearly 
less  than  half  the  height  H  of  the  pipe. 


We  shall  observe  more  symmetry  in  the  boundary  of  the  flow 
for  the  first  approximation.  The  lower  free  streamline  (L^)  is  found 
by  using  co^  of  case  1  above.  Hence  the  expression  is  ^ 


iw^x’ ) 
e  _ 

s/x"*"  (xf  -1  ) 


dx* 


(5.32) 


Under  the  change  of  variable 


x'  =  1/y1  we  obtain  ^ 

iCJ1  ( i/y* ) 

7= - 

s/77  (y’-l) 


-iw1(y' ) 

e 


dy 1 


) 


(5.33) 


1 

7t 


sfy7  (y'-l) 


-  73  - 


which  is  the  expression  for  the  upper  free  streamline  (L)  ,  providing 
the  exponent  in  the  integrand  be  changed  to  +  iu)^(y'  )  . 

This  means  that  the  upper  and  lower  free  streamlines  as  well 
as  the  surfaces  (s)  and  (S^)  are  symmetric  with  respect  to  a  line 
half  way  between  the  surfaces  (S^)  and  (s)  and  parallel  to  them. 
The  real  part  of  the  function  represents  the  slope  of  the  free 

streamlines  of  xT  >  a  or  0  <  x'  <  (3  .  In  table  four  is  recorded 
this  slope  for  various  values  of  r)  . 


Table  4 


T)  f1(n)  T)  fx(Tl) 


.05 

+ 

.003 

20.0 

-  .003 

.10 

+ 

.004 

10.0 

-  .004 

.20 

+ 

.005 

5.0 

• 

0 

0 

vn 

.50 

+ 

.011 

2.0 

-  .011 

.67 

+ 

.018 

1.5 

-  .018 

.91 

+ 

.075 

1.1 

-  .075 

•  95 

+ 

.181 

1.05 

-  .181 

.99 

+ 

.823 

1.01 

-  .823 

As  we  can  see  the  slope  becomes  zero  very  quickly  which  means  the  stream¬ 
lines  become  horizontal  very  quickly. 

The  approximations  of  order  greater  than  one  are  obtained  in 
a  manner  previously  outlined,  ^(z')  is  given  by  an  expression  similar 
to  (5.13)*  Hence  )  is  ^ 


■ 


-  ?4  - 


ii 

JC 


r sin  fi(Ti) 


-381(ti) 


Log 


_ in _ 

W  ^  (l-T,) 


(5-3^) 


where  D  *  [0,»)  -  [a,p]  ,  £1(n)  and  g,  ( rj)  are  obtained  from  case  1 

above  and  N(t),  zt)  is  defined  by  ^ 


(q-P)(zt-n)  +  (2n-a-p)  f(z,-|3)(a-z,)^+i(2zV-a-p)f(n-a)(il-p)}2  .  (5.35) 

(q-p)(z*  -rj)  +  (2r]-q-(B)  {(z‘  -p)  (q-z*  )  }^-i(2z 1  -q-0 )  { ( Tj-q)  ( q-0 )  }2 


This  must  be  interpreted  as  in  example  one.  Hence  | N ( rj, z * )  |  >  1 


by 


1 


The  lower  free  streamline  for  the  second  approximation  is  given 


1 

Jt 


iw1(x» ) 
e _ 

n/x*"  (x'-l) 


X  e 


(5.56) 


where  co^(x’)  is  obtained  from  case  1  above  and  £V,(x’  )  is  obtained 
from  (5.3^)* 


Now  the  real  part  of  (^(x*  )  will  represent  the  change  in 
the  slope  of  the  lower  free  streamline  for  the  first  approximation.  In 
the  interval  [0,0]  sin  f^  is  negative  while  in  the  interval  [q,00) 
it  is  positive.  Consequently  the  real  part  of  ^(x*)  for  x'  >  q 
will  be  negative. 


This  is  in  agreement  with  our  intuitive  argument  in  which 
gravity  acts  in  such  a  way  as  to  lower  the  lower  free  streamline  from 
the  first  approximation. 

In  table  5  is  recorded  a  typical  hand  computation  of  F^x1) 
for  x*  >  q  .  This  is  intended  only  as  an  indication  of  the  procedure 


in  evaluating  the  flow  boundaries.  One  finds  the  greatest  contribution 
of  the  integral  comes  firom  the  interval  [.4,3]  U  [a, 2]  since  outside 
of  this  sin  f^  is  almost  zero.  No  attempt  will  be  made  to  continue 
computation  on  the  second  approximation  or  to  proceed  to  the  third 
appoximation  without  the  use  of  a  computor  which  clearly  is  necessary. 


Table  5 


x?  =  15.0 


I 

n 

2 

3 

Log  |N(t],X»  )  | 

4 

Sin  f1(T]) 

5 

sT T)  (1-T]) 

6 

dl 

.02 

.04 

5.502 

-  .002 

.139 

-  .003 

.05 

.02 

5.475 

-  .003 

.212 

-  .002 

.07 

.02 

5.453 

-  .003 

.246 

-  .002 

.10 

.04 

5.422 

-  .004 

.285 

-  .003 

.20 

.  16 

5.311 

-  .005 

.358 

-  .012 

.40 

.22 

5.037 

-  .008 

.380 

-  .023 

.60 

.20 

4.644 

-  .014 

,310 

-  .042 

.80 

.15 

3.965 

-  .032 

.179 

-  .106 

.90 

.08 

3.280 

-  .067 

.095 

-  .185 

.95 

•  03 

2.591 

-  .137 

.049 

-  .217 

.97 

.02 

2.072 

-  .187 

.030 

-  .262 

.99 

.02 

.802 

-  .723 

.010 

-1.161 

1.01 

.005 

.800 

.723 

-  .010 

-  .292 

1.03 

.02 

2.070 

.24 

-  .030 

-  .320 

1.05 

.0425 

2.595 

.140 

-  .051 

-  .302 

1.10 

.125 

3.295 

.074 

-  .105 

-  .291 

1.30 

.20 

4.407 

.026 

-  .342 

-  .067 

1.50 

.30 

4.952 

.018 

-  .614 

-  .044 

2.00 

.80 

5 .664 

.011 

-1.414 

-  .035 

3.00 

1.00 

6.051 

.007 

-3.460 

-  .012 

4.00 

1.00 

6.929 

.006 

-6.000 

-  .007 

"6.00 

3.00 

7-641 

.005 

-12.25 

-  .009 

10.00 

5.00 

8.817 

.004 

-28.85 

-  .006 

14.00 

2.50 

11.51 

.003 

-48.64 

-  .002 

16.00 

2.50 

11.60 

.003 

-60.00 

-  .001 

25.00 

15.00 

7.54 

.002 

-120.00 

-  .002 

2  dl  =  -3.408 


sin  f  (t|) 

dl  =  Atj  Log  |N(rj,15)'|  - - - - 

V  T)  (  1  “  T) ) 


F0(15)  =  *  e"3Sl 


Z  dl 


. 00002  » 


For  more  results  see  Appendix  3. 
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CHAPTER  VI 

CONCLUSIONS 

Under  the  assumption  imposed  upon  the  flow  in  chapter  three 
we  have  investigated  the  flow  under  gravity  of  a  fluid  flowing  over  a 
bed  (S)  and  falling  freely  for  some  distance.  We  have  established 
the  uniform  convergence  of  the  sequence  of  the  successive  approximations 
to  a  function  which  satisfies  the  boundary  conditions.  We  have  out¬ 
lined  the  procedure  of  applying  these  results  to  concrete  flow  falls 
by  reference  to  two  examples. 

We  have  assumed  explicitly  that  an  equipotential  line  occurs  at 
a  point  downstream  of  the  point  of  overflow,  which  is  such  that  the 
variation  of  the  magnitude  of  the  velocity  along  it  is  zero.  In  actual 
fact  this  occurs  only  an  infinite  distance  downstream  but  it  has  been 
necessary  to  terminate  the  downstream  flow  domain  by  such  an  equi¬ 
potential  in  order  to  eliminate  a  theoretical  difficulty  which  other¬ 
wise  occurs.  This  assumption  may  be  satisfied  along  any  equipotential 
downstream  by  placing  sources  and  sinks  at  various  positions  outside 
the  downstream  flow  domain  but  this  affects  the  whole  flow  and  the 
origional  problem  is  changed.  The  effect  of  such  an  assumption  is  such 
that  the  nearer  the  equipotential  line  is  to  the  point  of  overflow 
the  more  serious  are  the  deviations  from  the  origional  problem.  Far 
downstream,  the  velocity  becomes  essentially  vertical,  so  from  the 
Bernoulli  equation  +  2gy  =  constant,  it  is  evident  that  the  deviation 

in  the  magnitude  of  the  velocity  along  every  equipotential  line  is 


■ 


. 
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very  small  and  the  effect  of  the  above  assumption  is  negligible. 

In  the  numerical  examples  that  have  been  examined  above,  the 
deviations  between  the  first  and  second  approximations  to  the  flow,  as 
are  evident  from  tables  3  and  4,  are  found  to  be  very  small.  This 
could  be  attributed  to  the  fact  that  the  equipotential  line  CC '  was 
too  close  to  the  point  of  overflow.  If  this  equipotential  line  is 
moved  a  greater  distance  downstream,  then  the  deviation  of  the  velocity 
along  this  equipotential  is  negligible  and  can  be  assumed  to  be  zero. 
This  will  make  some  difference  in  the  numerical  results  and  the  above 
discussion  points  to  the  fact  that  the  results  could  be  more  favourable. 
Due  to  the  small  size  of  the  results  we  have  obtained,  and  the  existence 
of  singularities  in  the  integrand  of  the  solutions,  accuracy  will 
necessitate  the  use  of  a  large  automatic  computor.  A  computor  study 
should  be  carried  out  before  definite  conclusion  can  be  drawn. 

The  pressure  can  be  determined  to  within  a  constant  at  any 
point  of  the  fluid  only  after  the  problem  is  completely  solved.  This 
constant  cannot  be  determined  because  the  origiona.1  pressure  p  ,  as 
it  occurred  in  the  Bernoulli  equation,  was  constant  along  the  free 
streamlines  and  does  not  appear  when  the  Bernoulli  equation  is 
differentiated.  In  assuming  that  an  equipotential  occurs  at  a  finite 
distance  downstream,  along  which  the  variation  of  the  velocity  is  zero, 
we  have,  it  appears,  assumed  that  outside  the  flow  domain  there 
occur  sources  or  sinks  which  make  this  possible.  It  is  not  possible, 
when  the  problem  is  solved,  to  determine  what  happens  outside  this 


flow  domain. 
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APPENDIX  I 


Continuity  of  Fn(l,s)  and  Gn(l,s) 


Define  $  (u)  =  sin  f  , (l,u)  e 
nv  '  n-lv  1  ' 


-5gn-i(1’u) 


R^(u)du  so  that 


Fn(i.s)  =  £ 

n  jt 


sin 

(u+s )/2 

sin 

CO 

1 

3 

du 


Since  R^(u)  is  zero  for  7  >  u  >  jt/2  the  integral  in 
Fn(l,s)  can  be  written  as  the  sum  of  two  integrals  whose  limits  of 
integration  extend  from  0  to  jt/2  and  from  7  to  it  .  These  two 
integrals  are  very  similar  since  the  type  of  singularity  and  order  of 
the  singularity  are  the  same.  We  shall  limit  ourselves  to  the  former 
of  these  two  integrals  since  continuity  of  the  latter  is  analogous  to 
that  of  the  former.  But  we  should  first  make  a  few  observations. 

Note  1,  Log  sin  u  is  Lipschitz  continuous  for  u  in  the  interval 
[a,b]  where  0  <  a  <  b  <  jt/2  . 


Note  2.  For  s  <  jt/2 


/  Log  sin  u  du |  =  |  /  (Log  u  + 
J  o  J  o 


,  sin  u  , 
Log  -  du 


u 


<  |  (Log  s-l)s |  +  |  / 

J  a 


_  sin  u 
Log  -  du 


u 


<  |s  |  {  |Log  s-1  |  +  |Log  | ) 

=  s  I  •  N  <  |  N 


’ 

' 


■ 


: 


*- 
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Note  3. 

s  i.n  A 

for  some  N  >  0  and  where  Log  — - —  is  the  maximum 

A 

value  of  Log  ^  for  0  <  u  <  jr/2  . 

The  integral  above  exists  since  all  of  the  singularities 

are  of  order  less  than  one. 

Note  4. 

$^(u)  is  continuous  in  the  interval  [0,d]  where  d  <  jt/2 

and  where  we  shall  assume  that  the  approximations  f^  ^(l,u) 

and  g  , (l.u)  are  bounded  and  continuous  on  the  circumference 
°n-lv  ’  ' 

of  the  unit  semi -circle. 

V1<s) 

The  following  cases  must  be  examined  in  order  to  prove  that 

is  continuous. 

Case  1. 

0  <  s  <  s 1  <  tf/2  y 

Case  2. 

jt/2  <  s  <  s'  <  ;t  ^ 

Case  3* 

s  jt/2“  v 

Case  4. 

s  -*n/2+? 

Case  5. 

s  0  ^ 

Case  6. 

V 

T 

(O 

We  shall  illustrate  the  procedure  by  considering  only  case  one. 

In  this  case  ^ 

o  <  ^  <  §  and  0  <  <  2  1  80 


. 

’ 

• 

... 


|Fn(l,.)-Fn(lf..)|  -  J  I  /  .n(“>{^  Iff Hf 


a)  |  i  $n(u){Log  si 

J  n 


U+S  T  .  U+S1  ->  J 

sm  -  Log  sin  — - —  J  du 


<  r  i$n(u)  jdu  x 

~~  J  n 


S-S1 I  <  € 


providing  js-s' |  <  6^  since  |®n(u) |  is  integrable  and  the  statement 
of  note  one. 


b) 


L 


u-s 


$>n(u){Log  sin  ,  2 


-  Log  sin  j~7”—j}  du 


/' 


$  (u)  {Log  sin  -  Log  sin 


s '  -u 


]  du 


s f 

+  J'  ®n(u)  [Log  sin  -  Log  sin  -p-)  du 


J '  ®n(u)  {Log  sin  ■  -  Log  sin  ~p—}  du 


=  | 2  P  [$  (s-2v)  -  $  (s’-2v)]  Log  sin  vdv 
J  n 


rt  8^ 

n8  “  4 

+  2  [®  (s+2v)-$>  (sf+2v)]  Log  sin  vdv 

Jo 


| -  j-  du 


. 
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it  s 

f ' * 

Jit  s* 

3  “  T 


$n(s+2v)  Log  sin  vdv 


Now  we  shall  deal  with  each  of  these  integrals  separately. 


<0 


f  {$n(s-2v) -®n(s ' -2v) )  Log  sin  vdv y 


<Cl*i  f  lL°8  sin  vl  dv  <  e  ) 


providing  js-a*  j  <  6^  and  by  notes  two  and  four  with  constant  >  0 
and  c  >  0  given. 


*) 


it  s 

3  “  T 


/  ($n(s+2v)«®n(s,‘j*2v) }  Log  sin  vdv  | 

J  o  n  7 


it 

3 


<  P  |®n(s+2v)-0n(sf+2v) |  X 


sin  v I  dv 


T 


<  c 


2  2 


J '  [Log  sin  v |  dv  <  e  ^ 


providing  Is-s1  I  <  6,  and  by  notes  two  and  four 

1  1  5 

and  e  >  0  given. 


constant  >  0 


®) 


p  e. 

/  $  (a**2v)  Log  sin  vdv| 

J «  / 


<  l  j§  (s’“2v)  Log  sin  vj  dv  <  Rjs-s’  j/2  <  e 
~  J  & 
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providing  js-s’  j  <  5^  for  c  >  0  given  and  using  the  fact  that  the 
integrand  is  bounded  in  this  interval. 


f) 


-2 


jt  sj_ 

pk  "  2 

Jt  S1 

H  "  IT 


0  (s'+2v)  Log  sin  vdv 


+  2 


n 


jt  £ 

Jt  s' 

8  "  4 


0  (s+2v)  Log  sin  vdv I 
n'  7  1  J 


it 


■<-L 


Jt  s' 

4  +  2 


v*»s  ^ 

<S>  (v)  Log  sin  — —  dv 
n  2 


X 


4  2 


+  s 


$  (v)  Log  sin  — r—  dv  1 
nv  '  2  1 


25.  +  11 
4  2 


.  /  n  rx  •  V-S  _  .  V-S 

®  (v;  [Log  sin  — Log  sin 


*  +  11 

pk  2  v_s 

+  /  $  (v)  Log  sin  -5—  dv! 

J  s'  jt  n  2 

S  "  T  +  4 


g) 


Jt 

,2 


25.  ,  £1 
4  2 


$  (v)  f  Log  sin 
n  L 


.  v-s 


Log  dv  I  % 

j 


< 


Jt 

,2 


i i  +  il 

4  2 


IVv) !  x 


l  X  K  dv  <  6  J 


dv 
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for  e  >  0  given  providing  |s-s*  |  <  &,_  and  using  note  one  and  the  fact 
that  |$  (v) |  is  integrable  in  the  interval  of  integration. 


h) 


Jt  s 1 

>4  +  2 


v-s 


.  0  (v)  Log  sin  — ~  dv I  <  |s-s*  I  Q  <  €  , 

s’  jt  nv  '  °  2  i=i  >  x  ? 


2  +  4 


providing  |s-sT  |  <  5^.  for  e  >  0  given  and  using  the  boundedness  of 
the  integrand  in  the  interval  of  integration. 


Then  for  e  >  0  there  exists  a  8  >  0  such  that  providing 
|s°sT  |  <  8  then  |Fn(l,s’)  -  F  (l,s)  |  <  €  for  0  <  s,s'  <  x/2  .  The 
remaining  cases  follow  in  the  same  manner  so  that  in  the  interval  [0,ot] 
Fn(l,s)  is  continuous,  A  similar  but  somewhat  simpler  proof  is 
necessary  to  prove  the  imaginary  part  0^(1, s)  is  continuous  in  this 
interval o 


APPENDIX  II 


b 

Suppose  we  wish  to  evaluate  the  integral  f  f(x)dx  and  f 

3 

cannot  be  integrated  directly.  Numerically  this  can  be  done  by  an 

approximation  procedure  known  as  the  method  of  tangents.  It  consists 

of  partitioning  the  interval  [a,b]  into  n  subintervals  [x^  ^,x_^] 

where  x  =  a  ,  x  =  b  and  Ax.  =  x.-x.  .  and  then  evaluating  the 
o  7  n  i  i  i“l 

expression 


(x) 


n 

i=l 


+ 


1 

2 


Ax .  )Ax , 

l  i 


Q 


If  f(x)  is  integrable  then  this  expression  can  be  made  to  approximate 

the  integral  as  closely  as  one  wishes  by  letting  maximum  {Ax.)  tend 

i  1 


to  zero. 
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APPENDIX  III 


Computation  of  Fg(x,o) 


Example  1 


X* 

F2(x‘) 

100. 

-  .061 

15. 

-  .037 

5- 

-  .020 

1.5 

-  .056 

1.05 

-  .103 

.95 

+  .085 

.05 

+  .041 

Example  2 


X*  - 

Fg(x') 

1.05 

-  .0009 

5.0 

-  .00045 

15.0 

-  .00027 
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